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,3 ■ ABSTRACT 

o 

(N ■ We present a quantitative model for the infrared emission from dust in the diffuse interstellar 

medium. The model consists of a mixture of amorphous silicate grains and carbonaceous grains, 
each with a wide size distribution ranging from molecules containing tens of atoms to large 
' grains > 1 micron in diameter. We assume that the carbonaceous grains have polycyclic aromatic 

hydrocarbon (PAH)-like properties at very small sizes, and graphitic properties for radii a > 50A. 
On the basis of recent laboratory studies and guided by astronomical observations, we propose 
^\ ' "astronomical" absorption cross sections for use in modeling neutral and ionized PAHs from the 

far ultraviolet to the far infrared. We also propose modifications to the far-infrared emissivity of 
"astronomical silicate" . 

We calculate energy distribution functions for small grains undergoing "temperature spikes" 
. caused by stochastic absorption of starlight photons, using realistic heat capacities and optical 

properties. Using a grain size distribution consistent with the observed interstellar extinction, we 
are able to reproduce the near-IR to submillimeter emission spectrum of the diffuse interstellar 
medium, including the PAH emission features at 3.3, 6.2, 7.7, 8.6, and 11.3/xm. The model is 
O I compared with the observed emission at high Galactic latitudes as well as in the Galactic plane, 

■ as measured by the COBE/DIRBE, COBE/FIRAS, IRTS/MIRS, and IRTS/NIRS instruments. 

^ ! The model has 60 x lO"'^ q (relative to H) locked up in PAHs, with 45 x IQ-^ of C in a 

component peaking at ~ 6A (Nq ~ 100 carbon atoms) to account for the PAH emission features, 
and with 15 x 10^^ of C in a component peaking at ~ 50A to account for the 60/xm flux. 
The total infrared emission is in excellent agreement with COBE/DIRBE observations at high 
' galactic latitudes, just as the albedo for our grain model is in accord with observations of the 

diffuse galactic light. The aromatic absorption features at 3.3/im and 6.2/im predicted by our 
dust model are consistent with observations. 

We calculate infrared emission spectra for our dust model heated by a range of starlight inten- 
sities, from 0.3 to 10** times the local interstellar radiation field, and we tabulate the intensities 
integrated over the SIRTF/IRAC and MIPS bands. We also provide dust opacities tabulated 
from the extreme ultraviolet to submillimeter wavelengths. 

Subject headings: dust, extinction — infrared: ISM: continuum — infrared: ISM: lines and bands 
— ISM: abundances — radiation mechanisms: thermal — ultraviolet: ISM 
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1. Introduction 



It has been 70 years since Trumpler's discovery of color excesses provided the first definitive proof 
of the existence of interstellar dust (Trumpler 1930); however, many aspects of the nature of interstellar 
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dust still remain unclear. Our knowledge of the composition of interstellar grains is based mainly on tlicnr 
spectroscopic absorption/emission features and observed elemental depletions. The most generally accepted 
view is that interstellar grains consist of amorphous silicates and some form of carbonaceous materials; the 
former inferred from the 9.7/im Si-0 stretching mode and 18/um 0-Si-O bending mode absorption features 
in interstellar regions as well as the fact that the cosmically abundant heavy elements such as Si, Fe, Mg 
are highly depleted; the latter mainly inferred from the strongest interstellar absorption feature - the 2200A 
hump - and the fact that silicates alone are not able to provide enough extinction. 

The dust sizes, inferred from the wavelength-dependent interstellar extinction and polarization curves, 
may be separated into two domains - 1) radii a > 0.02/im, which includes the "classical" grains (with a > 
0.1/xm), which arc primarily responsible for the extinction, polarization and scattering at visible wavelengths; 
and 2) the "very small grain" component (with a < 0.02/zm) which contribute importantly to the extinction 
in the vacuum- ultraviolet. While the size distribution for the "classical grain" component is relatively 
well constrained by fitting the observed interstellar extinction curve for an assumed dust composition, our 
knowledge of the size distribution for the "very small grain" component is much poorer, due to the fact that, 
for A > 0.1/um, these very small grains are in the Rayleigh limit and their extinction cross sections per unit 
volume are independent of size, so that the observed UV/far-UV extinction curve only constrains the total 
volume (mass) of this grain component. 

Forty-five years ago, Piatt (1956) proposed that very small grains or large molecules with radii < lOA 
may be present in interstellar space. Donn (1968) further proposed that polycyclic aromatic hydrocarbon- 
like "Piatt particles", may be responsible for the UV interstellar extinction. However, the inability to 
determine the detailed optical properties of "Piatt particles" and the limited observational information on 
the UV/far-UV interstellar extinction forestalled further development of these ideas (Greenberg 1960). 

Since the 1980s, an important new window on the "very small grain component" has been opened by 
infrared (IR) observations. The near-IR continuum emission of reflection nebulae (Sellgren, Werner, & Din- 
erstein 1983) and the 12 and 25/zm "cirrus" emission detected by the Infrared Astronomical Satellite (IRAS) 
(Boulanger & Perault 1988) explicitly indicated the presence of a very small interstellar dust component 
since large grains (with radii ~ 0.1/im) heated by diffuse starlight emit negligibly at such short wavelengths, 
whereas very small grains (with radii < 0.01/im) can be transiently heated to very high temperatures (> 1000 
K depending on grain size, composition, and photon energy) as first pointed out by Greenberg (1968). Sub- 
sequent measurements by the Diffuse Infrared Background Experiment (DIRBE) instrument on the Cosmic 
Background Explorer (COBE) satellite confirmed this and detected additional broadband emission at 3.5 
and A.dfi-ai (Arendt et al. 1998). 

More recently, spectrometers aboard the Infrared Telescope in Space (IRTS) (Onaka et al. 1996; Tanaka 
et al. 1996) and the Infrared Space Observatory (ISO) (Mattila et al. 1996) have shown that the difi'use 
interstellar medium radiates strongly in emission features at 3.3, 6.2, 7.7, 8.6, and 11.3/Ltm. These emission 
features, first seen in the spectrum of the planetary nebulae NGC 7027 and BD+30°3639 (Gillett, Forrest, & 
Merrill 1973), have been observed in a wide range of astronomical environments including planetary nebulae, 
protoplanetary nebulae, reflection nebulae, HII regions, circumstellar envelopes, and external galaxies (see 
Tielens et al. 1999 for a review). Often referred to as "unidentified infrared" (UIR) bands, these emission 
features are now usually attributed to polycyclic aromatic hydrocarbons (PAHs) which are vibrationally 
excited upon absorption of a single UV/visible photon (Lcger & Puget 1984; AUamandola, Tielens, & 
Barker 1985) although other carriers have also been proposed such as hydrogenated amorphous carbon 
(Duley & Williams 1981; Borghesi, Bussoletti, & Colangeli 1987; Jones, Duley, & Williams 1990), quenched 
carbonaceous composite (Sakata et al. 1990), coal (Papoular et al. 1993), fullerenes (Webster 1993), and 
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interstellar nano diamonds with sp^ surface atoms reconstructed to sp^ hybridization (Jones & d'Hendecourt 
2000). The emission mechanism proposed for the UIR bands - UV excitation of gas-phase PAHs followed 
by internal conversion and IR fluorescence - is supported by laboratory measurements of the IR emission 
spectra of gas-phase PAH molecules (Cherchneff & Barker 1989; Brenner & Barker 1989; Kurtz 1992; Cook 
et al. 1998) and by theoretical investigations of the heating and cooling processes of PAHs in interstellar 
space (AUamandola, Tielens, & Barker 1989; Barker & Cherchneff 1989; d'Hendecourt et al. 1989; Draine & 
Li 2001). 

By far the most extensively modelled dust property is the interstellar extinction curve. Existing grain 
models for the diffuse interstellar medium (see Witt 2000 for a recent review) are mainly based on an analysis 
of extinction (Mathis, Rumpl, & Nordsieck 1977; Greenberg 1978; Hong & Greenberg 1980; Draine & Lee 
1984; Duley, Jones, & Williams 1989; Mathis & Whiffen 1989; Kim, Martin, & Hendry 1994; Mathis 1996; 
Li & Greenberg 1997; Zubko 1999; Weingartner & Draine 2001a). The near-IR (l-S/xm), mid-IR (5-12/im) 
emission spectrum along with the far-IR (>12/im) continuum emission of the diffuse Galactic medium yields 
further insights into the composition and physical nature of interstellar dust; in particular, the PAH emission 
features allow us to place new constraints on the size distribution of the very small dust component. 

To model this IR emission, we require realistic calculations of the excitation and deexcitation rates of 
small grains subject to stochastic heating by absorption of starlight, and cooling by spontaneous emission 
of infrared photons. We require detailed knowledge of the optical and thermal properties of interstellar dust 
materials to derive the probability distribution for the vibrational energy of a grain of specified size and 
composition. 

Attempts to model the IR emission of interstellar dust have been made by various previous workers. 
Following the initial detection of 60 and lOO/xm cirrus emission (Low et al. 1984), Draine & Anderson (1985) 
calculated the IR emission from a graphite/silicate grain model with grains as small as 3 A and argued that 
the 60 and lOO/im emission could be accounted for. When further processing of the IRAS data revealed 
stronger-than-expected 12 and 25/zm emission from interstellar clouds (Boulanger, Baud, & van Albada 
1985), Weiland et al. (1986) showed that this emission could be explained if very large numbers of 3-loA 
grains were present. A step forward was taken by Desert, Boulanger, & Puget (1990), Siebenmorgen & 
Kriigel (1992), Schutte, Tielens, & AUamandola (1993), and Dwek et al. (1997) by including PAHs as an 
essential grain component. Early studies were limited to the IRAS observation in four broad photometric 
bands, but Dwek et al. (1997) were able to use DIRBE and FIRAS data. No calculations have been carried 
out to ascertain whether both the PAH features and the far-IR continuum can be quantitatively reproduced. 

There has been considerable progress in both experimental measurements and quantum chemical cal- 
culations of the optical properties of PAHs (AUamandola, Hudgins, & Sandford 1999a; Langhoff 1996; and 
references therein) . There is also an improved understanding of the heat capacities of dust candidate materi- 
als (Draine & Li 2001a, hereafter DLOl) and the stochastic heating of very small grains (Barker & Cherchneff 
1989; d'Hendecourt et al. 1989; DLOl) the interstellar dust size distributions (Weingartner & Draine 2001a, 
hereafter WDOla), and the grain charging processes (Weingartner & Draine 2001b). In the present paper we 
make use of these advances to model the full emission spectrum, from near-IR to submillimeter, of dust in the 
diffuse interstellar medium (ISM). The model is compared to IRTS, DIRBE, and FIRAS observations. We 
consider dust in three distinct regions: high-latitude regions (|6| > 25°; hereafter HGL); and two regions in 
the Galactic plane - the "MIRS" region (44° < / < 44°40', -0°40' < 6 < 0°) for which the 4.5-11. 7Aim spec- 
trum was obtained by the mid-IR spectrometer on board IRTS; and the "NIRS" region (47°30' < / < 48°, 
|6| < 15 ) for which the 2.8-3.9/im spectrum was obtained by the near-IR spectrometer on board IRTS. 
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In the present work we consider a dust model which is a natural extension of the original graphite-silicate 
model for interstellar dust (Mathis, Rumpl, & Nordsieck 1977; Draine & Lee 1984, hereafter DL84). Our 
dust model consists of two grain types - amorphous silicate grains and carbonaceous grains - both with 
size distributions ranging from ultrasmall grains/large molecules a few angstroms in size (in the molecular 
domain) to large grains > 1/im in diamctcr'^. Wc assume that the carbonaceous grains have graphitic 
properties at large sizes and PAH-like properties at very small sizes. For PAH-like grains, we distinguish 
between the optical properties of neutral PAHs and PAH ions. The current work assumes spherical grains 
and therefore does not yet address issues related to polarization. 

This paper is organized as follows. We first explore in §2 the optical properties of interstellar dust 

materials, with emphasis on the absorption cross sections of both ionized and neutral PAHs from far-UV 
to far-IR. We then briefly summarize the heat capacities of PAHs, graphite and silicates in §3. In §4 we 
describe the radiation field adopted for the diffuse ISM. We calculate in §5 the equilibrium temperatures for 
large grains (> 200A). In §6 wc describe the methods developed in DLOl to calculate the energy distribution 

functions for very small (a < 200A) grains heated by starlight. Selected results arc presented for illustration, 
and wc demonstrate the importance of stochastic heating for grains as large as a = 200 A. 

The adopted dust size distribution functions are summarized in §7. The PAH charging computation is 
described in §8, and §9 details the IR emission modelling. 

In §10 we discuss our choices of model parameters, and compare our model with observations of infrared 
emission from high galactic latitudes as well as two galactic plane regions. We find our model to be in 
very good agreement with observations from the near-IR (including the PAH emission features) to the 
submillimeter. We compute emission spectra for our "standard" dust heated by starlight with intensities 
ranging from 0.3 to 10^ times the local value, and we tabulate the predicted intensities averaged over the 
SIRTF IRAC and MIPS photometric bands. The total infrared emission from our grain model is seen to be 
in excellent agreement with DIRBE observations, just as the albedo for our grain model is in accord with 
observations of the diffuse galactic light. 

In §11 we show the extinction curve for our model, from the far-UV to submillimeter wavelengths, as 
well as tabulated values of the dust opacity. We discuss the elemental depletion, the aromatic absorption 
features, and the variations of the relative abundance of PAH in various regions. We also summarize the 
upper limits for the abundance of ultrasmall silicate grains (both amorphous and crystalline) estimated in 
Li & Draine (2001). Finally, we summarize our conclusions in §12. 

2. Optical Properties of Interstellar Dust Materials 

To calculate the energy (temperature) distribution functions of very small grains (undergoing "tem- 
perature fluctuations"), or the "equilibrium temperatures" of larger "classical" grains, and the resulting IR 
emission spectrum, an accurate knowledge of the absorption and emission properties of interstellar dust 
materials over a broad range of wavelengths is required. 



^To avoid confusion, we emphasize here that "very small grains", "large molecules", and "ultrasmall grains" are synonymous. 
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2.1. Silicate Grains 

We use Mie theory (see, e.g., Bohren & Huffman 1983) to calculate the absorption and scattering cross 
sections for amorphous silicate spheres. We adopt the dielectric functions of DL84 with two modifications: 
1) the sudden steep rise of the imaginary part in the far-UV has been smoothed (see WDOla); and 2) the 
dielectric function has been modified at A > 250/Ltm to better match the average high Galactic latitude dust 
emission spectrum measured by FIRAS (Wright et al. 1991; Reach et al. 1995; Finkbeiner, Davis, & Schlegel 
1999) - for A > 250/im we take the imaginary part of the silicate dielectric function to be 



e2(A) = e°i^(A) x 

where /3=0.4, and e^*^ is the imaginary part of the DL84 silicate dielectric function. Using the Kramers- 
Kronig relation (see DL84), we have recomputed the real part of the silicate dielectric function from the 
newly modified 62. The modification to 62 (A) and the absorption cross sections Cabs (A) for silicates is very 
slight except at very long wavelengths: there is no change at all for A < 250/zm, and for 250 < A < 1100/zm 
the new ei is within ±12% of . Only for A > 4.8 mm is the silicate emissivity altered by more than a 
factor of 2. 

We stress that we do not intend to suggest that the dielectric function for "astronomical silicate" 
has been determined with the precision that might be suggested by the modest adjustments represented 
by eq. (1) - i.e., we do not believe that €2 is known to an accuracy of ±12% for A < 1100/um. Our 
objective is simply to demonstrate that with an entirely reasonable dielectric function, our carbon-silicate 
grain model can accurately reproduce the available observations of the far-infrared emission spectrum. The 
original DL84 astro-silicate had Cabs A^^ ° for A > lOG/im. The above modification has Cabs « A~^'^ 
for 800 A < 1cm, similar to the A~^-^ frequency dependence observed by Agladze et al. (1996) over 
700-2500/(xm for amorphous 2MgO-Si02 grains at T = 20K [although note that at A = lOOO/xm we have a 
silicate opacity k — Cabs/(47r/9a^/3) = 0.33cm^g~^, whereas the material studied by Agladze et al. (1996) 
had K = 1.25 cm^ g~^ at this frequency]. 
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Q ln(A/250^m)xln(A/850;im) 
' P ln(850/250) 



(A/850/im)'' 
(10V850)'', 



250 < A < 850;tim, 
850 < A < lO^/xm, 
A > lO^/xm. 



(1) 



2.2. Carbonaceous Grains 

We assume a continuous distribution for the optical properties of interstellar carbonaceous grains with 
graphite-like properties for large sizes (> 50A) and PAH-like properties in the small size limit (< 20 A). The 
absorption cross sections CfJ^g(a,A) of graphite grains are calculated from Afie theory using the "1/3 2/3 
approximation" (Draine & Malhotra 1993) and the graphite dielectric functions of DL84. The absorption 
cross sections of PAHs will be discussed below (§2.3). Finally, we make a transition from PAH optical 
properties to graphite properties at a = 50A by taking the absorption cross sections of carbonaceous grains 
to be 

Ca^bf (a, A) = CPAHCr,^,«(a, A) + (1 - ePAH)Cf;;:(a, A) (2) 
CpAH(a) = (1 - ggra) X min [1, (aj/a)^] , = 50A, ggra = 0.01 (3) 
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where C^^, Cf[^g arc, respectively, the absorption cross sections of PAHs and graphite grains of radius^ 
a at wavelength A; is the grain radius where the transition from PAH properties to graphite properties 
begins; < ggra <C 1 so that even for small carbonaceous grains we introduce a small amount of "continuum" 
absorption (see §2.3 below); ^pah, the PAH "weight", drops from 1 — qg^a to as a increases from to 
infinity. Our choice of the functional form of ^pah is quite arbitrary, but the resultant IR spectrum is 
insensitive to the detailed behavior of ^pah provided ^pah = (1 — Qgra) for a < 30 A and ^pah ^ 0.5 for 
a > lOOA. For example, the much flatter function, ^pAH(a) = (1 ~ Qgva) x min {l, 1/ [l + (a/aj — 1)^] }, 
results in a model spectrum (not shown here) which is nearly indistinguishable from that of eq.(3). 



2.3. Absorption Cross Sections of PAHs 

In the present model, a PAH molecule is characterized by the number of carbon atoms (A^c), the 
hydrogen to carbon ratio (H/C), and charge state (neutral or charged). The H/C ratio is determined by the 
aromatic structure. Since it seems likely that astronomical PAHs are relatively pericondensed, as favoured by 

stability considerations (van der Zwct & Allamandola 1985; Omont 1986; AUamandola et al. 1989; Jochims 
et al. 1994) and the fact that large compact PAHs provide a better match to interstellar PAH band positions 
(Langhoff 1996), we assume 

r 0.5, A^c < 25, 

H/C = <^ 0.5/v'A^c/25, 25 < A'c < 100, (4) 
[ 0.25, A^c > 100, 

which for A^c < 100 approximates those of compact, symmetric PAHs (series 1, 3 in Stein & Brown 1991). 
But we have also considered a higher H/C ratio which is appropriate for PAHs with more open, uneven 
structures (see §10.6). 

Following previous workers (Leger, d'Hendecourt, & Dcfourneau 1989a; Desert et al. 1990; Schutte 
et al. 1993), we construct the absorption cross sections for both neutral and ionized PAHs from far-UV 
to far-IR based on available laboratory measurements and guided by the astronomical observations (see 
Appendix A for details). The resulting cross sections are characterized by a set of Drude profiles: the cr-cr* 
transition with a peak at Aq"^ ~ 14/im^i; the tt-tt* transition with Aq"^ ~ 4.6/im^i; the C-H stretching mode 
(Ao = 3.3/xm); two C-C stretching modes (Aq = 6.2,7.7/xm); the C-H in-plane bending mode (Aq = 8.6/um); 
three C-H out-of-plane bending modes (Aq = 11.3, 11.9, 12. 7/(im); and a few weak features probably due 
to C-C bending modes (Aq = 16.4, 18.3, 21.2, 23. ly^m). In addition to these discrete features, there is also 
far-UV, near-UV/ visible and far-lR continuum absorption. 

The absorption cross section per C atom C^^/Nc is taken to be 

CtllW/Nc X > 17.25; (5) 

(126.0 - 6.4943a;) X 10^1® cmVC, 15 < a; < 17.25; (6) 

S'i(A) (-3.0-M.35a;) X 10-i^cmVC, 10 < a; < 15; (7) 

(66.302 - 24.367a; + 2.950x^ - 0.10573;^) x lO'^^ cmVC, 7.7 < a; < 10; (8) 

^We adopt a mass density ai 2.24 g cm^'' for graphite. The "radius" a of a PAH containing Nq C atoms is defined to be 

1/3 S 

the radius of a sphere with the carbon density of graphite containing the same number of C atoms, i.e., a = 1.286Nq A. 



A^c 
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= S'2(A)+ 1.8687 + 0.1905a; + 0.4175 (a;- 5.9) V 

0.04370 (x - 5.9f 1 x 10"^^ cmVC, 5.9 < a; < 7.7; (9) 

= 52(A) + (1.8687 + 0.1905x) x 10"^^ cm^/C, 3.3 < x < 5.9; (10) 

14 

= 34.58 X 10-18-3.431/a; ^ cutoff (A, Ac) cmVC + ^ S'j (A), X < 3.3; (11) 

i=3 

where x = (A//Ltm)~^. The ultraviolet a-a* and n-n* transitions (j=l,2), C-H stretching (j=3), C-C stretch- 
ing (i=4,5), C-H in-plane bending (i=6), C-H out-of-plane bending (j=7,8,9), C-C bending (j=10,ll, 12,13) 
modes, and far-IR continuum (j = 14) are represented by Drude profiles, with the absorption cross section 
per carbon atom contributed by feature j given by 



2 7jAj(7int.j 



sjw ^ ''T:t:. , ,2 . (12) 



where Tabic 1 lists values of the central wavelength Xj, the broadening parameter 7j, and the integrated 
absorption strength^ 

aintj = J S,{X)dX-' = ls,iX,)^,X-\ (13) 

The function cutoff(A, Ac) and the cutoff wavelength Ac, determined by the PAH size (the number of fused 
benzcnoid rings), arc defined in cqs. (A2,A3). 

Figure 1 plots the absorption spectrum of PAHs in the UV/far-UV. Also plotted are the experimental 
spectra of coronene (C24H12) (Leger et al. 1989b), PAH mixtures (Leger et al. 1989b; Joblin, Leger, & Martin 

1992), and the cross section estimated for graphite (DL84). Figure 2 displays the absorption cross sections 
of neutral and ionized PAHs for a wide wavelength range with emphasis on the C-H, C-C IR vibrational 
modes (see following and §2.4, §A.3, §A.4 for details). 

We have included enhancement factors £^6.2, -£^7.7, and £^8.6 in the strengths of the 6.2, 7.7, and 8.6/xm 
bands^. For £'6.2 = -B7.7 = -Eg. 6 = 1 we recover the laboratory values, but we are unable to satisfactorily 
reproduce astronomical spectra unless we take £^6.2, £^7.7! £■8.6 > 1- For the 7.7 and 11.3/im features, 
laboratory data (see §A.3) indicate a relative band strength crint.7.7/cint.ii.3 ~ 1-3 for neutral PAHs, and 
ft! 12 for ionized PAHs (assuming H/C = 1/3). Let P{Xj) = J IxdX be the total power integrated over 
the profile of feature j. Values of P(7.7/im)/P(11.3/im) as large as 11 are observed in some objects, such 
as compact HII regions (Roelfsema et al. 1996). DLOl have modelled the emission expected from PAHs in 
various radiation fields, and find that for £^7.7=2, fully ionized PAHs with iVc <1G0 in a region of very intense 
UV radiation (up to lO^x the Habing [1968] UV intensity) achieve P{7.7iJTa)/P{11.3inji) « 21. However, 
since the observed spectra include emission from larger PAHs with smaller values of P(7.7/im)/P(11.3/im), 
it is necessary to take £^7.7 w 2 to reproduce the strength of the strongest observed 7.7/um features. 

In addition, a wide variety of objects show a P(6.2/xm)/P(7.7/im) ratio exceeding the maximum pre- 
dicted from the laboratory PAH feature strengths. For example, P(6.2/xm)/P(7.7/im) « 0.3 - 0.6 in a 



^The cross section integrated over wavelength J Sj{\)d\ = A|(Tintj- 

^Absolute band strengths for some PAH species vary considerably among different laboratory groups and between experi- 
mental results and theoretical calculations. For example, the integrated cross sections of naphthalene, anthracene, tetracene, 
phenanthrene, benzanthracene, pyrene, and coronene cations can differ by a factor of 1.5-20 between different experimental 
studies, or between experimental and theoretical studies (see Langhoff 1996). 
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"With £'6.2=1, Er. 7=1 and Es.6 = l one obtains the laboratory measured band strengths. The observed IR spectrum of the 
diffuse ISM is best fitted by our dust model with -E6.2=3, i?7.7=2 and £8.6=2- 

'For the C-H out-of-plane bending modes, the j factor arises from the assumption that PAHs have equal numbers of H in solo, 
duet, and trio adjacent CH units (sec §A.3). 
''Band width from Moutou et al. (2000). 
''Band widths from Moutou et al. (1996). 

''To reproduce the weak X > 14/im continuum absorption measured by Moutou et al. (1996). 
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Fig. 1. — The ultraviolet and far ultraviolet absorption cross section per C atom of PAHs. 
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Fig. 2. — Absorption cross section per C atom of neutral and ionized PAHs (with Eq,2=3, ^7.7=2, Eg. 6=2). Ionization 
enhances the C-C stretching modes (6.2, 7.7/im) and C-H in-plane bending mode (8.6/xm), and weakens the C-H stretching 
mode (3.3/im). Ionization shifts the visual absorption edge (which is also size-dependent) redward. The absorption cross section 
per C atom of small graphite spheres is shown for comparison. 
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recent study of 820 mid-IR spectra of the diffuse interstellar medium (Chan et al. 2001), and « 0.54 for 
an average of planetary nebulae, reflection nebulae, and HII regions (Cohen et al. 1989). Thermal emission 
has a strict upper limit P(6.2/;im)/P(7.7;um) < (7.7/6. 2)2(7i„t,6.2/o-i„t, 7.7 « QAAEQ,2/E^,^ (see Table 1), and 
detailed modelling of stochastic heating finds P(6.2/im)/P(7.7/im) < 0.35£'6.2/-E'7.7 (see Figure 16 in DLOl). 
Therefore, we take Eq,2 = 1-5-B7.7 « 3. 

For most objects, the strength of the 8.6^m C-H in-plane bending mode is consistent with the strength 
of the 7.7/im C-C bending mode, so we take iSg.e = -E7.7 w 2. 



2.4. Continuum Opacity of Polycyclic Aromatic Hydrocarbons 

Even when a < aj, eqs. (2,3) cause Cabs to include a fraction ^gra of the opacity per C atom of graphite, 
resulting in "continuum" opacity extending into the infrared. This is required to account for the near- 
IR continuum emission at 1 5/im underlying the 3.3/im C-H stretching feature detected in a number of 
astronomical objects, including the Orion Bar (Joblin et al. 1996a; Bregman et al. 1989; Geballe et al. 1989), 
various reflection nebulae (Sellgren, Werner, & AUamandola 1996), the Red Rectangle (Geballe et al. 1989), 
normal galaxies (Helou et al. 2000), and the starburst galaxies M82 and NGC 253 (Sturm et al. 2000). This 
continuum emission may arise from the electronic fluorescence or a quasi-contimmm of overlapping overtone 
and combination bands from vibrational fluorescence in PAHs (Leger & Puget 1984; AUamandola et al. 
1985, 1989; Sellgren et al. 1996). Electronic intraband transitions can also contribute an optical-infrared 
continuum due to the closing of the band gap as the PAH size increases and/or upon ionization of the PAH 
(AUamandola et al. 1989; Tielens et al. 1999). 

We represent this continuum by the parameter ^g^a in eq. (3). High resolution spectroscopic observations 
of NGC 1333 (Joblin et al. 1996a), the Orion Bar (Geballe et al. 1989; Jobhn et al. 1996a), and the Red 
Rectangle (Geballe et al. 1989) resolving the 3.3/im feature with FWHMpa 0.04/um (same as that adopted 
here) show a feature peak-to-continuum ratio as high as ^ 30. If we require a 3.3/im peak-to-continuum 
ratio of 30, then PAHs with A^c ^ 50 (small enough to be heated to T > 1500 K by a UV photon - see 
DLOl) must have 

_ 1 2 X 3.3/im X aint,3.3Mm _ n ni /tA\ 



where we have taken C^^g(3.3/tm) « 3x 10~^°iVc cm^, H/C« 1/3, and have taken (Tint.s.s^im to be the average 
of the values in Table 1 for neutral and ionized PAHs. The resulting PAH infrared cross sections are shown 
in Figure 2, for both neutral and ionized PAHs. 

Finally, there is the question of the value of a^. The observed PAH emission features are at A < 25/im, 
and this will be contributed primarily by particles with a < 40A which can be heated to > 100 K by a 
photon with hc/X < 13.6 eV. Therefore, the observed PAH emission features only require > 25 A. We 
will take w 50A so that the 6.2/im feature in absorption is in agreement with observations, as discussed 
in §12, but this value - as well as the form of the transition from PAH- like to graphitic optical properties 
(§2.2) - is poorly determined, as the calculated emission spectrum is relatively insensitive to variations in 
in the range 25 - 100 A. IR spectra calculated with = 25A or = 75A (not shown here) are almost 
identical to that for = 50A. 
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3. Enthalpies 

As discussed in DLOl, the vibrational energy of a PAH molecule containing Nc C and (H/C) x Nq H 
atoms can be approximated as 

3 

EpAu{T) = ^Nc J2 e^p{ti^/kT) _ i + (^c - 2) [kQopf2{T/eop) + 2keMT/&^p)] ; (15) 

the summation index j = 1 — 3 runs over the C-H out-of-plane bending (a;i/27rc=886 cm~^), in-plane 

bending (a;2/27rc=1161 cm"-'^), and stretching (ti)3/27rc=3030 cm"'^) modes; the term containing Nc — 2 is 
a two-dimensional Debye model with Debye temperatures Qop ~ 863 K for the out-of-plane C-C modes, 
@ip w 2504 K for the in-plane C-C modes; k is the Boltzman constant; T is the "vibrational temperature"; 
and 

y"dy 



fn{x) = - f 



(16) 



/o cxp(y/x)-l" 
For graphitic grains, we use eq. (15) with H/C = (see DLOl). 
For amorphous silicate grains, DLOl recommend 

E,^{T) = (iVatom - 2) [2fce2/2(r/e2) + fce3/3(r/e3)] , (i7) 

where A^atom is the number of atoms contained in the silicate grain (we assume a mass density 3.5g cm^"^ 
for MgFeSi04 silicate material); 62 ~ 500 K, 63 « 1500 K are Debye temperatures for the two- and 
three-dimensional Debye models, respectively. 



4. Radiation Field 

We adopt the solar neighbourhood interstellar radiation field (ISRF) (Mathis, Mezger, & Panagia 1983; 
hereafter MMP) for the diffuse ISM 

C 4 



ux = XMMP <! + V Wi—Bx{T,) \ + — Sa(2.9K) (18) 



c c 

i=2 ) 



where Xmmp describes the enhancement of the starlight component (relative to the solar neighborhood)^; 
u^^® is the ultraviolet component; {W2, W3, W4) = (10"", 10-^^ 4x lO'^^^ and (Ta, T3, T4) = (7500, 4000, 3000) 
K are dilution factors and corresponding black-body temperatures. The last term is the cosmic background 
radiation. 



5. Thermal Equilibrium Temperatures of Large Grains 

Large grains maintain a nearly constant temperature T determined by balancing absorption and emis- 



sion, 



pOO POO 

/ C^bs{a,X)cuxdX= / Cabs(a,A)47rBA(T)rfA , (19) 
Jo Jo 



^Notc that Xmmp = 1 corresponds to x = 1-23, where x is the ratio of the intensity relative to the Habing (1968) estimate 
at lOOOA, or Go = 1.14, where Go is the 6-13.6 eV energy density relative to the 6-13.6 eV energy density in the Habing (1968) 
radiation field. 
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where CabsCojA) is the absorption cross section for a grain with size a at wavelength A, c is the speed of 
hght, B\{T) is the Planck function at temperature T, and u\ is the energy density of the radiation field (see 
§4). In Figure 3 we display these "equilibrium" temperatures for graphitic and silicate grains as a function 
of size in environments with various UV intensities. 



6. Transient Heating of Very Small Grains 

Since Grcenberg (1968) first noted that very small grains heated by starlight would experience temper- 
ature spikes, there have been a number of theoretical studies of this process. DLOl have recently carried out 
a detailed study in which they compared the energy distribution functions and resulting emission spectra 
calculated "exactly" (the exact-statisticaF method) versus two approximations: the "thermal- discrete" ap- 
proximation and the "thermal-continuous" approximation. The exact-statistical method involves calculation 
of the vibrational density of states. In contrast, both the thermal-discrete method and the thermal-continuous 
method assume that at any given time, the vibrationally excited grain has a "vibrational temperature" , and 
the deexcitation process is assumed to correspond to a thermal emission spectrum. In the thermal- discrete 
method both the excitation and deexcitation processes are treated as discrete events; however, in the thermal- 
continuous method, the deexcitation process is approximated as continuous rather than discrete. 

Although the thermal methods are not exact, DLOl have shown that the overall IR emission spectrum 
calculated using the thermal-discrete and the thermal- continuous methods is very close to that of the exact 
treatment (see also Barker & Chcrchneff 1989; d'Hendecourt ct al. 1989; AUamandola ct al. 1989; Schutte 
et al. 1993). Since the exact-statistical method involves heavy computations of the grain vibrational mode 
spectrum and density of states, we adopt the thermal- discrete method. Figures 4 and 5 show the energy 
distribution functions found for PAHs and silicate grains with radii a = 5, 10, 25, 50, 75, 100, 150, 200, 300A 
illuminated by the average ISRF (xmmp = 1)- Very small grains (a < 100 A) have a very broad P{E), 
and the smallest grains (a < 30 A) have an appreciable probability Pq of being found in the vibrational 
ground state = 0. As the grain size increases, P{E) becomes narrower, so that for xmmp = 1 it can be 
approximated by a delta function for a > 250 A. However, for radii as large as a = 200 A, both carbonaceous 
and silicate grains have energy distribution functions which are broad enough that the emission spectrum 
deviates noticeably from the emission spectrum for grains at a single temperature T, as shown in Figure 
6. For accurate computation of infrared emission spectra it is therefore important to properly calculate the 
energy distribution function P{E), including grain sizes which are large enough that the average thermal 
energy content exceeds a few eV. 



7. Grain Size Distributions 

Following Weingartner & Draine (2001a), we adopt size distribution functions for silicate grains and 
carbonaceous grains of the following form: 

1 dnsii Csii f a ^ 



riH da a Vat 



F{a; (3t,sii,at) x G(a; at,sii, ac,sii) , (20) 



rf^carb ^ j^^^^ Ccarb / « \ F(a; /3carb, Ot) X G(a; flt.carb, ^ccarb) , (21) 



uh da 

^ V^t.carb 



^^'*'^'"*^ = M/(l-/3aM), /3<0 ' 
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Fig. 3. — Equilibrium temperatures for graphite (dotted lines) and silicate grains (solid lines) in environments with various 
starlight intensities. 
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E/hc (cm-i) 

Fig. 4. — The energy distribution functions for charged carbonaceous grains (a = 5 [C6oH2o]i 10 [C4soHi2o]i 25 [C720oHi80o]i 
50, 75, 100, 150, 200, 250, 300A) for Xmmp = 1- The discontinuity in the 5, 10, and 25A curves is due to the change of the 
estimate for grain vibrational "temperature" at the 20th vibrational mode (see DLOl). For 5, 10, and 25A a dot indicates the 
first excited state, and Pq is the probability of being in the ground state. 
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Fig. 5. — As in Figure 4, but for silicate grains. 
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X (/xm) 

Fig. 6. — Infrared emission spectra for small carbonaceous grains of various sizes heated by starlight with xmmp = 1, 
calculated using the full energy distribution function P{E) (solid lines); also shown (broken lines) are spectra computed for 
grains at the "equilibrium" temperature T. Transient heating effects lead to significantly more short wavelength emission for 
a < 200 A. 
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, X / 1> 3.5A<a<at , . 

G{a;at,ac) = < c \/ i3\ ^ ' (23) 

[ cxp |-[(a- at)/acj''| , a > at 

where Uh is the hydrogen number density, and D{a) is an additional population of very small (a < 100 A) 
carbonaceous grains which is introduced in order for the grain model to reproduce the observed infrared 
emission. The grains in the distribution D{a) are all in the Rayleigh limit at the wavelengths A > 1000 A 
where interstellar extinction has been measured, and therefore extinction measurements constrain only be, 
the number of carbon atoms per total H atoms in the population D(a). For the PAH ultraviolet absorption 
cross section discussed in §2.3, WDOla showed that values of be < 60ppm are allowed for the average diffuse 
ISM extinction law {Ry = 3.1). For a given choice of 6c within the allowed range, the parameters at,sih 
flc^sii, ttsii, /3sii, Csii, at,carb, Oc.carb, ^carb, /3carb, C'carb Can be found in WDOla. 

The population of very small grains required to explain the IR emission of the cirrus and translucent 
molecular clouds is larger than could be obtained from simple extrapolation of the MRN power law {dn/da oc 
a-^-^) to very small sizes (Weiland et al. 1986; Desert et al. 1990; Dwek et al. 1997; Verter et al. 2000), and 
therefore D{a) must include such very small grains. The observed infrared emission requires a substantial 
population of grains in the 5 10 A size range so that single-photon heating can raise them to high enough 
temperatures to emit in the 3 — 12/xm range (see DLOl). A log-normal distribution was used by Draine & 
Lazarian (1998) to represent the size distribution of ultrasmall carbonaceous grains producing both 3 — 12/im 
infrared vibrational emission and microwave rotational emission, and we adopt such a distribution here. 

We will see below (see Figure 8) that the a > 100 A grains in the present grain model, with the 
adopted optical properties of graphite and silicate grains, produce less emission than observed near 60/im 
for a > 100 A grains, as the steady-state grain temperatures T are slightly too cool. This could be because 
the grain properties may be in error - perhaps Cabs is too low in the optical and ultraviolet, or too large 
in the far-IR, or perhaps the starlight radiation field has been imderestimated. However, we find that 
the grain model gives excellent agreement with the total radiated power (suggesting that the product of 
starlight intensity and optical- UV absorption cross sections are about right), and excellent agreement with 
the A > lOO/xm emission (suggesting that the far-IR emission cross sections are about right). 

As previously suggested by Draine & Anderson (1985), the shortfall at 60/im could be made up by 
emission from a population of very small grains which single-photon heating can raise to ^ 40 K, so that 
much of the absorbed energy will be radiated near 60/Lim. The optimal grain size for this purpose is a ~ 50 A. 
We use a second log- normal component to add a population of a « 50 A grains to D{a): 



D{a) 




ln(a/aoi'"' ^ 



cr 



a > 3.5 A (24) 



p _ 3 exp(-4.5g'^) mcbcA ^ , . 

' (2^)3/2 al^pa {l + erf[3a/y2 + ln(ao./3.5A)/a^/2]} ' ^'"^'^ ^ 

where mc is the mass of a C atom; p = 2.24 gcm"^ jg the mass density of graphite; aoi determines the 
location of the peak, bc,i is the amount of C atoms relative to H consumed by the i-th log-normal com- 
ponent. The width of each of the log-normal components is determined by a. The lower cutoff is set at 
o=3.5A (corresponding to Nq w 20 for PAHs) since smaller grains are photolytically unstable (Omont 1986; 
AUamandola et al. 1989; Guhathakurta & Draine 1989). Component i of D{a) has a^dn/dhia peaking at 

Ctpeak.i = Ooiexp(3cr^). 

As we will show below, we obtain a good fit to the observed interstellar emission spectrum with bci = 
45ppm, aoi = 3.5 A, bc2 = 15ppm, ao2 = 30 A, and a = 0.4. The resulting size distribution, with a^dn/dlna 
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having three separate peaks at ~ 6 A, ~ 50 A, and 3000 A (see Figure 2 of WDOla) seems somewhat artificial, 
but this may be attributable in part to the procedure used for fitting the UV extinction. 



8. PAH Charging 

Due to their low ionization potential (6-7 eV), PAH molecules can be ionized through photoelectric 
emission (AUamandola et al. 1985; van der Zwet & AUamandola 1985; Lepp & Dalgarno 1988; Bakes & 
Tielens 1994). On the other hand, PAHs may also acquire charge through collisions with electrons and 
ions (Draine & Sutin 1987). The grain charging process has been recently reanalyzed by Weingartner & 
Draine (2001b); we use their rates for photoelectric emission and electron capture to study PAH charging in 
three idealized interstellar environments characterizing the diffuse ISM: cold neutral medium (CNM), warm 
neutral medium (WNM), and warm ionized medium (WIM). The parameters adopted in the PAH charge 
distribution computation are displayed in Table 2. 

We estimate (/'ion (a) - the probability of finding a PAH molecule of radius a in a non-zero charge state 
- for CNM, WNM, and WIM respectively. In Figure 5 we present the PAH "ionization fraction" ^ion(a) as 
a function of size for CNM, WNM, and WIM conditions, and in addition the average 0ion(a) for the mass 
fractions listed in Table 2. 



9. Modeling the IR Emission Spectrum 

With the absorption cross sections Cabs (a, A) (§2), the UV radiation field strength xmmp (§4), the grain 
size distribution (§7), the energy distribution functions P{E) for very small grains (§6), the equilibrium 
temperature T for large grains (§5), and the ionization fraction (pion{a) of PAHs (§8), we can now calculate 
the IR emission spectrum from a mixture of silicate and carbonaceous grains. The IR emissivity per H 
nucleon, for dust illuminated by starlight of strength x, is 

Exix) = I r'^7da—^CtL{a,X)Bx{fsn[a,x]) 
[J250A da 

+ da^^ClZ{a,X)B,if,,4a,xd 
J250A '^H da 

+ r^da^^ r dE(^^) CtL{a,X)B,{T,^[EM) 

+ r"^da—^^ r di;Uo„ f^) Cl^:^\a,X)B^{T,,,^^[E,a]) 



Table 2: Idealized phases for the diffuse ISM . 



Item 


CNM 


WNM 


WIM 


Mass fraction 


0.43 


0.43 


0.14 


nH(cm-3) 


30 


0.4 


0.1 


T(K) 


100 


6000 


8000 


rie/riH 


0.0015 


0.1 


0.99 



20 



20 50 102 IQ3 C IQ4 




a (A) 

Fig. 7. — The "ionization fraction" 0ion — the probabiUty of finding a PAH molecule in a non-zero charge state - as a function 
of size for three idealized phases of the diffuse ISM: CNM (dotted), WNM (short-dashed), and WIM (long-dashed). Solid line 
shows the averaged ionization fraction weighted by the mass fractions of Table 2. The slope discontinuity at a = 50A is due to 
the transition from PAH optical properties to graphite properties above = 50 A. The upper axis label give the number of 
carbon atoms A^c = 0.468(a/ A)^. 
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+ (1 



i^cart 



(a, X)Bx{T^arb'' [E,a]) 



87r/ic/A5 



(26) 



where C|}Jg, Cf^^, C^^* , and C^§s^° are the absorption cross sections for silicate grains, graphitic grains, 
and neutral and ionized carbonaceous grains, respectively; Tsu, Tgra are the equilibrium temperatures of 
silicate and graphite grains, respectively; dP{E) is the probability that the grain vibrational energy will 
be in [E,E + dE]; and Tx{E,a) is the effective vibrational temperature for a grain with composition x, 
vibrational energy E, and radius a (see DLOl). The first term on the right represents the contribution of big 
silicate grains (a > 250A; hereafter "Bgii"); the second term represents the big carbonaceous grain component 
(a > 250A; hereafter "Bcarb")i third term represents the very small silicate component (3.5A< a < 250A; 
hereafter "Ssii"); the fourth and fifth terms represent the charged and neutral very small carbonaceous grains 
(including PAHs) for 3.5A< a < 250A; hereafter "Scarb"- The ux/{S'!Thc/X^) factor is a (usually negligible) 
correction for stimulated emission. In regions where dust is optically thin to its own thermal emission, the 
resulting IR intensity is simply 

(A)dust = Exix) X Nn (27) 

where A^h is hydrogen column density. Dust self-absorption should be taken into account in regions with 
large Nu- 

(T \ _ p r^W 1-exp [-AHSabs(A)] 

^abs(.AJ 

where Sabs('^) is the total absorption^ cross section per H nucleon for the dust model. 



10. Model Emission Spectra and Observations 
10.1. Overview of Model Parameters 



If we choose a value for a (we take a = 0.4), then the grain model has a total of 18 adjustable parameters: 



1. gTgra, determining the PAH near-IR continuum (1-5/im) opacity; we take ^gra = 0.01 (see §2.2, §2.4); 

2. a^, determining the transition from PAH optical properties to graphite properties; we take = 50A 
(see §2.2); 

/3carb, Ccarb, characterizing the size distributions and 
abundances of silicate and carbonaceous grains (§7); 

4. aoi, bc,i, ao2, &c,2 (= be — bcs), characterizing the size distributions and abundances of the two 
log-normal components of very small carbonaceous grains (§7); 

5. Xmmp, the starlight intensity (§4); 

6. Ah, the total column density of material. 



In addition, the spectrum depends on the ionized fraction ^ion(a) for PAHs, determined by XMMp/^^e {ne is 
the electron number density) and the gas temperature T. 



If the radiation field in the emitting region is approximately isotropic, scattering does not alter the intensity. 
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For the UV properties of PAHs adopted here, WDOl show that the observed interstehar extinction puts 
an upper limit on 6c = 6c, i + bc,2 ^ 60 ppm for the average extinction in diffuse regions (WDOl). For a 
given choice of 6c < 60 ppm, ten of the parameters (at,sii, ac.sii, ccsii, Psii, Csh, at,carb, ac,carb, "carb, Pcarh, 
Ccarb) are determined by fitting the interstellar extinction curve (see WDOl a). 

The strength of the 5-12/im IRTS spectrum relative to the total far infrared emission in the MIRS 
region requires a population of PAHs with 6c, i ~ 45 ppm with a mass distribution peaking near ~ 6 A. Thus 
Ooi ~ 6 Aexp(— 3ct^); adopting a = 0.4 gives aoi ~ 3.5 A. 

For the high galactic latitude dust we take xmmp = 1, since the starlight heating this nearby, optically- 
thin diffuse material should be very close to the local ISRF estimate of Mathis, Mczgcr, & Panagia (1983). 

The dust at high galactic latitudes has substantially more 60/zm emission than would be produced 
by graphite and silicate grains at their "thermal equilibrium" temperatures T. Within the context of the 
graphite-silicate model, we can account for this extra emission with a population of very small carbonaceous 
grains which undergo single-photon heating to T w 40 K, which corresponds to a « 50 A grains. We therefore 
add 6c. 2 ~ 15ppm in a population of grains with ao2 ~ 50 Aexp(— 3(t^) « 30 A. 

This completes specification of the properties of our "standard" grain model, with parameters summa- 
rized in Table 3. To compute the emission, it remains only to specify the intensity of the starlight heating 
the grains, and the total column density A^h- 



10.2. Observational Constraints 

We will compare our model results with observational data for 3 regions: the average emission per 
H observed by DIRBE (Arendt et al. 1998) and FIRAS (Finkbeiner et al. 1999) for high-Galactic-latitude 
regions; the DIRBE flux and the 4.5-11.7/im IRTS spectrum (Onaka et al. 1996) for the "MIRS" region; 
and the DIRBE flux and the 2.8-3.9/im IRTS spectrum (Tanaka ct al. 1996) for the "NIRS" region. For 
the MIRS and NIRS regions the DIRBE fluxes were obtained from the DIRBE data set, obtained from the 
National Space Science Data Center (NSSDC). In Table 4 we list the DIRBE fluxes for these three regions. 



Table 3: Parameters for standard dust model for Ry = 3.1 



PAHs 



carbonaceous*^ 



sihcate** 



6c,i=45ppm, aoi=3.5A (ct=0.4) 
&C,2=15ppm, ao2=30A (c7=0.4) 
aj=50A, Qgra = O.Ol 
-£'6.2 = 3, -B7.7 = Eg Q = 2^ 
H/C from eq. (4) 



Ot,carb = 0.0107/um 
flc.carb = 0.428/im 
CKcarb = -1.54 
/Jcarb = -0.165 



flt.sii = 0.164/im 
flc.sii = 0.100/Ltm 
Qsii = -2.21 
/?sii = 0.30 



Ccarb = 9.99 X 10-12 c'sil = 1.00 X 10-13 



"Taken from Wcingartner & Draine (2001a). 

*See Table 1. 
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10.3. Model Spectra: High-Galactic-Latitude Emission 

The high galactic latitude regions are optically thin to starlight, so the starlight illuminating the grains 
is approximately uniform, with intensity Xmmp ~ 1- We adopt our standard grain model, with parameters 
from Table 3. 

Arendt et al. (1998) have determined A/^/AiYn, which wc compare with E\{x) from cq. (27). In Figure 
8 we show XIx/Nn for our dust model for xmmp = 1- The model is in excellent agreement with the observed 
dust emission per H nucleon for A > lOO/xm. At 60/im the model emission is ~ 18% low, which is probably 
within the uncertainties in the DIRBE fluxes^ and the uncertainties in subtraction of zodiacal emission in 
high galactic latitude regions (Arendt ct al. 1998). There is a significant discrepancy between model and 
observation at 12/xm, where the model flux is only 55% of the observed flux. In view of the difficulties in 
correction for zodiacal emission, which is relatively strong at 12/im, the shortfall at 12/im may not be real. 
However, agreement at 12/xm could be improved if the PAH ionization fraction (/'ion(a) were reduced for 
6 <a< 10 A, since in this size range the strong 7.7/xm feature in ionized PAHs radiates away energy at the 
expense of the 11.3, 11.9, and 12.7/um C-H out-of-plane bending modes. Unfortunately, there are as yet no 
observational data regarding the strength of the 6.2, 7.7, and 8.6/im emission from high galactic latitudes. 

The model emission also appears to fall short at S/zm, where the model flux is only 40% of the observed 
value. At these shorter wavelengths correction for zodiacal emission is increasingly difficult. However, if the 
reported AIx/ANn is correct it would require additional opacity at this wavelength for the ultrasmall grain 
component. Such additional opacity could, for example, be contributed by carbon chain molecules with a 
carbon abundance of < Ippm (Allamandola et al. 1999b). 

The far-IR-submillimeter emission from dust at high galactic latitudes has been measured by FIRAS 
(Wright et al. 1991; Reach et al. 1995; Finkbciner et al. 1999). In Figure 9 we show the average FIRAS 
spectrum from Finkbeiner et al. (1999) together with the far-IR emission calculated for our grain model, 
for grains illuminated by xmmp = 1- Even using the original DL84 dielectric functions for graphite and 
silicate, the model is in very good quantitative agreement with the FIRAS measurements for A < llOO/xm. 
The agreement is improved if the silicate cmissivity is modified for A > 250/im, as described in §2.1. The 
resulting emission spectrum is then in excellent agreement with the observations. Finkbeiner et al. (1999) 
had argued that the FIRAS spectrum was indicative of two types of dust grains, one of which was quite 



Dwek et al. [1997] have estimated that the DIRBE fluxes have 1-cr uncertainties of ±20%. 



Table 4: DIRBE results for high galactic latitudes and NIRS and MIRS regions . 



Item 


2.2^m 


3.5/xm 


4.9/im 


12/im 


25/im 


60/im 


100/im 


140yum 


240/im 


Aeffectivo (/^m) 


2.22 


3.53 


4.88 


12.3 


20.8 


56.0 


97.7 


148 


248 


FWHM (/im) 


0.37 


0.97 


0.67 


8.1 


8.9 


27.3 


31.0 


35 


101 


HGL"* 




0.97 


1.11 


7.16 


3.57 


5.30 


18.6 


22.5 


10.1 


MIRS region'' 


6.65 


3.47 


1.68 


3.79 


2.10 


4.85 


11.5 


19.4 


8.20 


NIRS region^ 


7.13 


3.47 


1.55 


3.01 


1.58 


3.63 


9.01 


14.9 


6.27 



''A(A/A)/AAfH for |6| > 25° (IQ-^e ergs"! sr-iR-i) (Arendt et al. 1998). 

^XIx (lO-^ergs-icm-^sr-i) for 44° < Z < 44°4o', -0°4o' < 6 < 0° (Hauser et al. 1998). 

"XIx (lO-^ergs-icm-^sr-i) for 47°30' < / < 48°, |6| < 15 (Hauser et al. 1998). 
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Fig. 8. — Comparison of the model to the observed emission from the diffuse ISM at high galactic latitudes > 25°). 
Curves labelled and B^arb show emission from "big" (a > 250 A) silicate and carbonaceous grains; curves labelled Sgii and 
Scarb show emission from "small" (a < 250 A) silicate and carbonaceous grains (including PAHs). Triangles show the model 
spectrum (solid curve) convolved with the DIRBE filters. Observational data are from DIRBE (diamonds; Arendt et al. 1998), 
and FIRAS (squares; Finkbeiner et al. 1999). 
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cold. However, we see here that the "graphite-sihcate" grain model, with the sihcate emission modified by 
no more than ±12% for A < 1100/zm, gives essentially perfect agreement with the FIRAS spectra. 



10.4. Model Spectra: Galactic Plane Regions 

The DIRBE short wavelength channels are dominated by starlight. We model the DIRBE spectra by 
taking 

^ 1 - exp [-iVHSabs(A)] 



-'ab! 



:(A) 



(29) 



where (/A)dust is obtained from eq. (28), the multipUcativc factor A determines the strength of the starUght, 
and the Wi and Ti are the same as for the ISRF (see §4). Equation (29) would be exact if the stellar density 
was proportional to the dust density and the starlight was isotropic at all points along the line of sight (so 
that scattering docs not alter the intensity). While not expected to be highly accurate, this approximation 
should give a reasonable approximation to the overall spectrum. 

We adopt our "standard" grain model (see Table 3). This assumes that the dust in the 6 kpc molecular 
ring (where most of the interstellar matter in the MIRS and NIRS fields is located) has the same size 
distribution as dust in the local diffuse ISM; this assumption is likely to be incorrect in detail, but is the 
best we can do at this time. 



10.4- 1- Single-x Models 

The grains in the MIRS and NIRS Galactic plane regions will be illuminated by a range of radiation 
fields, since each 0.5° x 0.5° field (52 pcx52 pc @ 6 kpc) will include diffuse clouds, molecular clouds and 
star-forming regions. The simplest approximation is to assume that the dust is uniformly distributed in the 
beam and heated by a uniform radiation field with intensity xmmp relative to the local ISRF. For a given 
Xmmp, the column density Nh is adjusted to match the DIRBE MOi^m intensity. The constant A in eq. 
(29) determining the strength of the stellar contribution has been adjusted so that the spectrum agrees with 
the DIRBE 3.5/im measurement. We find that the Xmmp=2 model provides a rough overall fit to the IRTS 
spectra and DIRBE photometry. 

In Figure 10 we show the emission spectrum calculated for the "standard" dust model, for grains heated 
by Xmmp = 2, for both the MIRS and NIRS regions. Note that the model spectra have been smoothed to 
the resolution of the NIRS and MIRS instruments on IRTS for comparison with the IRTS spectra. Also 
note that the column densities A^h are very large, corresponding to Ay = 23 and 17 mag for the MIRS and 
NIRS regions, respectively. These column densities are large enough so that self- absorption is significant 
at A < 15/Ltm. In particular, the 3.3, 6.2, and 11-12/im features are noticeably attenuated by the dust 
absorption. 

This simple model provides fair agreement with the DIRBE observations of these two fields. The largest 
discrepancy is at lOO/im, where the observed flux is only 70% of the model value. 

The model is in reasonably good agreement with the IRTS 5-12/xm spectrum in the MIRS region. The 
model is ~40% below the 11.3/im peak observed by IRTS, and ~32% below the DIRBE 12^m measurement, 
but note that 1) the claimed absolute calibration accuracy of the MIRS instrument is only ±30% (Onaka et 
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Fig. 10. — Infrared emission from dust plus starlight (see text) for two regions in the Galactic plane: (a) the MIRS region 
(44° < I < 44°40', -0°40' < b < 0°), and (b) the NIRS region (47°30' < I < 48°, \b\ < 15'). For both regions we use our 
standard dust parameters, with the PAH ionization fraction <t'ion{o-) taken to be the weighted mean for the local ISM (43% 
CNM, 43% WNM, 14% WIM). The starlight intensity heating the dust has been taken to be Xmmp = 2. The solid curve shows 
the overall model spectrum, degraded to the IRTS resolution; triangles show the model spectrum convolved with the DIRBE 
filters. DIRBE observations are shown as diamonds. For the MIRS field we show the IRTS MIRS 5-12;im spectrum (thin solid 
line). For the NIRS field we show the IRTS NIRS 2.8-3.9^m spectrum (thin solid line, also shown as cross-dotted curve in 
inset). Curves labelled Bsji and B^arb show emission from "big" (a > 250 A) silicate and carbonaceous grains; curves labelled 
Ssii and Scarb show emission from "small" (a < 250 A) silicate and carbonaceous grains (including PAHs). 



28 



al. 1996), 2) the DIRBE photometric accuracy is estimated to be ±20% (Dwek et al. 1997). 

The model is also in reasonably good agreement with the IRTS NIRS 2.8 3.9/im spectrum in the NIRS 
region. Recall that the stellar contribution has been adjusted to match the DIRBE 3.5/im photometry; with 
this choice of stellar continuum, the model falls <15% below the NIRS spectrophotometry for A < 3.7/xm. 
This agreement is satisfactory, since the difference of the overall slope of the 2.8-3.9/zm NIRS spectrum and 
the slope indicated by the DIRBE 2.2, 3.5, and 4.9/im photometry suggests that the uncertainties in the 
NIRS spectrophotometry may exceed the claimed ±5% absolute calibration uncertainty (Tanaka et al. 1996). 
The model prediction for the strength of the 3.3/tim emission feature is in good agreement with the NIRS 
photometry (see the inset in Figure 10b). 

In the context of the single-x model, dust models with either weaker or stronger radiation fields are less 
successful than the xmmp=2 model in reproducing the observed spectra. The xmmp=1 model does not emit 
sufficiently at 5 12/im due to self- absorption by the very large column density of dust; on the other hand, 
the xmmp=3 model, although providing a reasonably good fit to the 5-12/L(m spectrum, emits too much (by 
> 75%) at lOOyum. 



The overprediction of the IGG/im flux in Figure 10 could be due in large part to the simplification of 
assuming that all of the dust is heated by a single radiation field. Consider instead a model where a fraction 
fj of the dust is heated by starlight with intensity Xj- If the regions of different Xj Si,ve randomly located on 
scales small compared to the beamsize, then we may approximate the different dust emissivities as uniformly 
mixed, so that the emergent intensity becomes 



The xmmp = 2 model of Figure 10 exceeds the observed 100/im emission by about 40%. This might indicate 
that the dust used to reproduce the 140/im emission is somewhat too hot. We can try to remedy this by 
using cold dust to provide the 140 and 240/im emission, plus additional warmer dust to provide sufficient 
60/im emission. Two examples of multi-x models for the MIRS region are shown here. In Figure 11a we 
show a multi-x model where the cooler dust is heated by xmmp = 0.7, and the warmer dust by xmmp = 100. 
This model reproduces the 60, 100, 140, and 240/im photometry, but the mid-IR absorption by the large 
amount of cool dust causes the 11-12/um emission to be too low by more than a factor of two. If the cooler 
dust is taken to be heated by xmmp = 1, we obtain Figure lib. This model now exceeds the observed lOO/xm 
flux by ^20%. The total dust column density has been reduced by ~30% compared to Figure 11a, but is 
still ~80% greater than in the single-x model of Figure 10a - there is sufficient absorption at 12/im that this 
model falls short of the observed flux by 30%. We have also carried out calculations for models combining 
three or more x values and found no significant improvement. We conclude that multi-x models offer no 
significant improvement over the simple single-x model of §10.4.1. 



10.4.2. Multi-x Models 




(30) 



10.5. Emission Spectra for Different Starlight Intensities 



In §10.3 we have shown models for optically-thin emission from our standard model heated by starlight 
with Xmmp = 1 for the high-latitude dust. To study the sensitivity to the assumed starlight intensity, in 
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Fig. 11. — Same as Figure 10a, but for multi-x models. Upper panel: a model where the far-IR emission is provided by 
dust with Xmmp = 0-7, with additional emission from dust heated by Xmmp = 100. While providing good agreement with the 
60-240/xm DIRBE fluxes, the model is a factor of two low in the 12/im region, partly due to the absorption contributed by the 
very large column density of Xmmp = 0.7 dust. Lower panel: a model where the far-IR emission is provided by dust heated by 
Xmmp = li with additional emission from dust with Xmmp = 100. There is less dust absorption, but the model still falls short 
in the 12/jm region. 
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Fig. 12. — As in Figure 8, but for dust illuminated by starlight intensities Xmmp = 0.5 (upper panel) or xmmp = 1-5 (lower 
panel) . 
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Figure 12 we show spectra for xmmp = 0.5 and 1.5. As expected from Figure 8, the xmmp = 0.5 model 
yields too Httle emission (by « 50%) over the entire spectral range; the Xmmp=1-5 model emits too much 
at A > 50/zm. 

Finally, in Figure 13 we show the emission from dust heated by radiation fields with xmmp = 0.3, 1, 3, 10, 100, 10"^, lO** 
(we assume all have the same ^ion as the local ISM). We see that the fraction of the infrared emission emerg- 
ing at A < 30/im is relatively independent of Xmmp (over this range), since this part of the spectrum is 
dominated by single-photon heating. This naturally explains why the observed PAH emission features show 
little spectral change in various regions with UV intensities differing by five orders of magnitude (e.g., see 
Boulanger 1999). However, as expected, the emission at A > 50/um shows a shift from longer to shorter 
wavelengths as xmmp is increased from 0.3 to 10^. 

The Space Infrared Telescope Facility (SIRTF) will be capable of sensitive imaging using the Infrared 
Array Cam,era, (IRAC) at 3.6, 4.5, 5.8, and 8.0/im, and using the Multihand Im,aging Photometer (MIPS) at 
24, 70, and 160/xm. In Table 5 we show the band-averaged intensities for our standard dust model heated 
by starlight radiation fields of various intensities. 

Table 5: I„/Nu (MJysr-VlO^^cm-^) Averaged over SIRTF Bands for Standard Dust Model Illuminated 
by Starlight Intensity xmmp- 
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IRAC 




IRAC 


IRAC 


MIPS 


MIPS 


MIPS 


Xmmp 


3.6/^m 


4.5/xm 




5.8/um 


8.0/xm 


24/xm 


70/xm 
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3.49x10-"* 
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0.188 


1.200 


3 
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0.0726 
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4.51 
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10^ 
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1690 
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10.6. 


PAH Feature Strengths 


and the H/C Ratio 





Our standard model assumes enhancements Eq,2,i.7,9,.& = (3,2,2) of the PAH band strengths. To 
demonstrate the need for such enhancement, in Figure 14a we plot the best-fitting model using the laboratory- 
measured PAH feature strengths [i.e. i?6.2,7.7,8.6 = (1, 1, !)]• Although the fit to the 11.3/im peak becomes 
a little better than our "standard" model (Figure 10a), the fit to the 6.2, 7.7, 8.6/im region is not as good 
as the standard model. We also note that the fit does not necessarily get improved by increasing the PAH 
feature strengths by a larger factor as expected from energy conservation consideration. 

Previous sections all assume a H/C ratio representing compact, symmetric PAHs (eq. [4], §2.3). For 
PAHs with open, uneven structures (e.g., elongated PAHs), the H/C ratio is expected to be higher. With 
this in mind, we consider a dust model with a relatively higher H/C ratio 



H/C = 



r 0.5, 

I 0.5/^/iVc760, 
[ 0.25, 



Nc < 60, 
60 < A^c < 240, 
Nc > 240. 



(31) 
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Fig. 13. — IR spectra (scaled by the starlight intensity) for our standard grain model for different illuminating radiation 
fields. The overall shape for A < 15/im is almost unchanged for regions with starlight intensities varying by about five orders of 
magnitude. Note we have assumed the same ionization fraction {</>ion) as in the local diffuse ISM. We also show the wavelength 
coverage of the 4 IRAC bands and 3 MIPS bands on SIRTF. 



33 



The corresponding spectrum is shown in Figure 14b. The C-H features (3.3, 8.6, 11.3/im) become sUghtly 
stronger, and the C-C features (6.2, 7.7/xm) become shghtly weaker (as expected from energy conservation), 
but the entire spectrum is nearly indistinguishable from the model with a lower H/C ratio (eq. [4]; Figure 
10a). 

The dust model described in this work assumes essentially full hydrogen coverage for PAHs. Dehydro- 
gcnation becomes significant only for the smallest PAHs. Upon absorption of an energetic photon, large 
PAHs will not have CH bond rupture since the absorbed energy will promptly be redistributed among many 
vibrational modes. Tielens et al. (1987) concluded that small PAHs (iVc < 25) would be partially dehydro- 
genated in regions with intense UV fields, while large PAHs (A^c ^ 25) would be completely hydrogenated 
in those regions. AUain, Leach & Sedlmayr (1996) reexamined this and reached similar conclusions - as 
shown in their Figure 10, PAHs with Ac > 30 in diffuse clouds are almost 100% hydrogenated. Therefore, 
it is reasonable to assume full hydrogen coverage. Figure 14b shows that the infrared emission spectrum is 
insensitive to modest changes in the H/C ratio (except for possible effects on the wavelengths and strengths 
of the C-H, C-C bands upon dehydrogenation [e.g., see Pauzat, Talbi, & EUinger 1997]). 

In regions very rich in H atoms, a PAH molecule may acquire excess H atoms attached to the peripheral 
C atoms and become "superhydrogenated" . The addition of extra H atoms to PAHs results in a decrease 
of the level of "aromaticity" and an increase of the level of "aliphaticity" by weakening the strength of the 
aromatic C-H stretching band and creating new aliphatic C-H stretching bands (Bernstein, Sandford, & 
Allamandola 1996). The 3.40, 3.46, 3.51, and 3.56/im subfeatures superposed on the shoulder of the 3.3/im 
feature seen in many objects may be explained by such "superhydrogenation" (Schutte et al. 1993; Bernstein 
et al. 1996). The very strong 8.6^m band seen in the compact HII region IRAS 18434-0242 (Roelfsema et al. 
1996) may be due to the extra H atoms bonded to internal C atoms, although the presence of non-compact 
PAH ions can also enhance the 8.6/im band (Roelfsema et al. 1996). Due to the lack of sufficient experimental 
data and quantum chemical calculations, we have not included possible "superhydrogenation" in our model. 

10.7. Total Infrared Emission and Optical-UV Albedos 

For Xmmp = 1 our model for the high galactic latitude emission (Figure 8) yields a total infrared 

intensity Att J E\d\ = 4.36 x 10"^'' ergs s^^ H^^ (the fractional contributions of Scarb, Bcarb: and silicate 
grains [Sgii+Bgii], are approximately 33%, 37%, and 30%, respectively). As seen from Figures 8 and 9, our 
model appears to be in excellent agreement with the overall power radiated by interstellar dust. 

By contrast, the model of Dwek et al. (1997) radiates ^37% more power per H nucleon (6.0x 10~^^ erg s~^ H~^), 
and the model of Desert et al. (1990) radiates ^59% more power per H nucleon (6.92 x 10~^^ergs~^ H~^). 
These differences arc due to the fact that 1) the model of Dwek et al. (1997) has a UV extinction which is 
higher than the average Galactic extinction curve (see their Figure 4) and probably has lower optical-UV 
albedos than ours; 2) the albedos predicted by the Desert et al. (1990) model are significantly lower than 
those of the current model. 

In Figure 15 we show the albedos calculated for the present model (also see Figure 15 in Weingartner 

& Draine [2001a]) together with the dust albedos estimated from observations of the diffuse Galactic light 
(Lillie & Witt 1976; Morgan, Nandy, & Thompson 1976; Hurwitz, Bowyer, & Martin 1991; Witt, Friedmann, 
& Sasseen 1997). The current model appears to be in good agreement with the observationally-determined 
albedos. For comparison. Figure 15 also shows the albedo for the grain model of Desert et al. (1990), 
which evidently has an absorptivity about 20-50% higher than both the current model and observations at 
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A < 4/xm where the bulk of the starlight energy is located. 

11. Discussion 

Figure 16 compares the extinction curve obtained by our dust model with 6c = 60ppm (6c, i = 45ppm, 
bc,2 = 15ppm) with the average interstellar extinction curve and also shows the contributions of the indi- 
vidual dust components. It can be seen that PAHs, the dominant contributor to the 2175A hump, have 
a relatively flat behavior for 5 < A^^ < 8/im^^ and thus do not dominate the far-UV extinction. This is 
consistent with Greenberg & Chlewicki (1983) in which the 2175A hump was shown to correlate poorly with 
the far-UV extinction. 

Our dust model implies a total dust mass per H nucleon mdust ~ 1-89 x 10~^^ g H~^, or a gas-to-dust 
ratio 1.4mH/TOdust = 124. Let Sext(A) = N^^ J{dn/da)Cext{0'^ ^)da be the total extinction cross section 
per H nucleon. The mass absorption coefficient for the dust is just Kabs('^) =^ (1 ^ £^lbedo)I]ext(A)/rndust = 
0.4 InlO (1 - albedo) [>1(A)/A^h] /mdust = 4.88 x 10^5(1 - albedo) [A(A)/7Vh] cm^ g-^. In Table 6 we tabulate 
Kabs at selected wavelengths. Note that these numbers are for dust in the average diffuse ISM, with Ry = 
Av/E(B — V)k, 3.1. In diffuse regions with Ry > 3.1 the extinction per H is reduced at A < 0.9/im due 
to changes in the grain size distribution (see WDOl), but the extinction does not appear to be affected 
at A > O.di^m (Cardelli, Clayton, & Mathis f989). In dark clouds ice mantle formation adds additional 
absorption features, most prominently the 3.1/im 0-H stretch in H2O, but for A > 10/um the opacities in 
Table 6 probably remain a good approximation even for molecular regions with densities up to 10^ cm~^. 

Our dust model requires a total of 254ppm of C and 48 ppm of Si in grains, both of which are significantly 
higher than the maximum available for the dust according to the recent work of Snow & Witt (1995), 
who suggest that the elemental abundances in the ISM are just 60%-70% of the solar value^. However, 
Weingartner & Draine (2001a) have argued that it is difficult to see how the observed interstellar extinction 
can be explained if interstellar abundances arc appreciably subsolar; they argue that the case for subsolar 
abundances is not compelling. Dwek (1997) has argued that the composite dust model (Mathis 1996) - 
which aims at making economical use of dust materials - is still short of carbon when PAHs are included to 
account for the mid-IR emission. 

Our model has an absorption feature at 6.2/Ltm due to the PAH C-C stretching mode (see Figure 16). The 
integrated optical depth of the 6.2/im feature relative to the silicate optical depth is (J AT(iA^"'^)6. 2/^^9.7 ~ 
2.9 cm~^. The aromatic 6.2/xm absorption feature has been detected in several objects including both local 
sources and Galactic center sources (Schutte et al. 1998; Chiar et al. 2000). On average, the observed 
integrated strength of the 6.2/im band per unit silicate band depth is about 2.2cm~^ with no significant 
difference between the Galactic sources and the local sources (Schutte et al. 1998). Our model is in good 
agreement with this determination. 

In addition to the 6.2/Lim feature, our model also predicts a weak narrow absorption feature at 3.3/Lim 
with J Ars.adA^^/Arg.T ~ 0.50 cm~^. This feature has been detected in the Galactic center source GCS 
3 with (J ArdA~^)3.3/AT9.7 ~ 0.34cm~^ (Chiar et al. 2000). It has also been seen in absorption in some 
heavily extincted molecular cloud sight lines (Sellgren et al. 1995; Brooke, Sellgren, & Geballe 1999). The 



'^Although the present model uses more silicon than other existing models, the predicted silicate optical depth AT9.7^m/-'VH ~ 
3.3 X 10~^^ cm^/H is not inconsistent with observations (Arg.rpm/iVH » 2.8 — 5.9 X 10~^^ cm^/H; see Draine 1989b, Mathis 
1998). 
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2 4 6 8 10 



A-1 (/xm-1) 

Fig. 15. — (a): Albedo calculated for the present dust model (solid; also see Figure 15 in Weingartner & Draine 2001a) and 
for the Desert et al. (1990) model (dotted; see their Figure 6) together with observational data for the diffuse Galactic light 
(filled circles — Lillie & Witt 1976; open squares - Morgan et al. 1976; open circles — Hurwitz et al. 1991; triangles - Witt et 
al. 1997). (b): The solar neighbourhood (xmmp = 1) interstellar radiation field (Mathis et al. 1983). The Desert et al. grain 
model is much more strongly absorbing at < 4fim^^ where the ISRF energy is concentrated. 
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0.1 1 10 100 1000 

A (/i-m) 

Fig. 16. — The calculated extinction for the present grain model (heavy solid line) and the observed average 
{Rv ~ 3.1) interstellar extinction curve (dotted line: Fitzpatrick 1999; open circles: Mathis 1990). Model 
results are the sum of 4 components: "Bsu" (silicate, a > 250A; long-dashed line); "Ssu" (silicate, 3.5 < a < 
250A; short-dashed line); "Bcarb" (carbonaceous, a > 250A; dot-dashed line); and "Scarb" (carbonaceous, 
3.5 < a < 250A, including PAHs; thin solid line). The upper-right inset illustrates the aromatic 6.2/im 
absorption feature and the 9.7/im silicate band. The narrow feature at 11.5/im is due to a lattice mode of 
crystalline graphite (Draine 1984) and would likely be smoothed out in an imperfect polycrystalline sample. 
The lower- left inset plots the extinction curve against (/xm"^). 
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Tabic 6. Extinction and absorption properties for dust in the diffuse ISM.^ 
NOTE: values of Sext and Kabs in ApJ Table 6 are too small by a factor of 1.086 . 



band 


A 


albedo 


g = (cos d) " 


^ext ^ 


d 

^abs 




ifim) 






/ 2 TT — 1 \ 

(cm'' H '^j 


/ 9 — 1 \ 

(cm^g 1) 




0.0100 


0.47 


1.00 


7.11 X 10—^^ 


2 no X 1 o4 




0.0488 


0.30 


0.86 


1 52 X 10—21 


5 66 X 104 




0.0912 


0.24 


0.73 




1.01 X 10'^ 




0.100 


0.27 


0.72 


2 28 X 10-21 


8 8^^ X 1 fl'^ 




0.1215 


0.32 


0.73 


1.73 X 10~^^ 


U . A _LU 




0.150 


0.38 


0.70 




4 56 X in4 




0.220 


0.42 


0.56 


1 65 X 10—21 


5 05 X in4 




0.300 


0.58 


0.57 


1 on V 1 n~2i 


9 9^ y HT^ 


u 


0.365 


0.62 


0.58 




1 fi8 X ID'^ 

X « yjKJ /N A.KJ 


B 


0.440 


0.65 


0.57 


6 73 X 10-22 


1 26 X in4 


Y 


0.550 


0.67 


0.54 


Ofi V 1 n~22 

O.UU A ±U 


8960 


R 


0.700 


0.66 


0.48 




6460 


J 


0.900 


0.63 


0.40 


2.47 X 10~^^ 


4820 


J 


1.22 


0.58 


0.29 


1 X in~^^ 


3440 


H 


1.63 


0.51 


0.21 


Q 76 X 1 0—23 


2510 




2.20 


0.43 


0.13 


'^ Q9 V 1 n~23 

iJ.t?Zi A ±U 


1780 


L 


3.45 


0.28 


0.0051 


2 57 X 10-23 


984 




3.60 


0.26 


-0.0037 


2 37 X 1 0—23 


926 




4.50 


0.18 


-0.037 


1.54 X 10—23 


670 


M 


4.80 


0.16 


-0.043 


1.37 X 10—23 


611 




5.80 


0.099 


-0.052 


1 02 X 10—23 


488 




8.00 


0.016 


-0.050 


1.54 X 10-23 


803 




9.70 


0.003 


-0.036 


4.39 X 10-23 


2320 


N 


10.6 


0.003 


-0.032 


3.36 X 10-23 


1770 




12.0 


0.003 


-0.029 


1.93 X 10-23 


1020 


Q 


21.0 


5 X 10"* 


-0.024 


1.30 X 10-23 


693 




25.0 


3 X 10-4 


-0.024 


9.29 X 10-24 


493 




60.0 


1 X 10-* 


< 0.01 


1.91 X 10-24 


101 




70.0 


1 X 10-" 


< 0.01 


1.38 X 10-24 


73.1 




100 


< 10-4 


< 0.01 


6.40 X 10-25 


34.0 




140 


< 10-4 


< 0.01 


3.31 X 10-25 


17.6 




160 


< 10-4 


< 0.01 


2.44 X 10-25 


13.0 




240 


< 10-4 


< 0.01 


1.06 X 10-25 


5.61 




450 


< 10-4 


< 0.01 


2.80 X 10-26 


1.49 




850 


< 10-4 


< 0.01 


8.83 X 10-2^ 


0.468 




1350 


< 10-4 


< 0.01 


4.05 X 10-2'^ 


0.215 




2000 


< 10-4 


< 0.01 


2.10 X 10-2'^ 


0.111 




3000 


< 10-4 


< 0.01 


1.05 X 10-2^ 


0.0559 


"An 


expanded 


version 


of this 


table is 


available a 



http : //www . astro . pr inceton . edu/~draine/dust /dustmlx . html 



''Average of cos d for scattered light, where d is the scattering angle. 

•^Sext is the dust extinction cross section per H nucleon. 

'^Kabs is the absorption opacity of the dust alone; the opacity of the gas-|-dust 
mixture is Kabs/125. To within ±10%, Kabs 3.17 X 105 (A//im)-2 (cm2 g-l) 
for 20< A<700/im, and Kabs " 3.89 x 104 (A/zim)-'-®® cm2 g-i for 700< A< 
104/im. 
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7.7, 8.6, and 11.3/im features are hidden by the much stronger 9.7/iin siUcate feature, and therefore will be 
difHcult to observe as absorption features. 

The aromatic absorption bands allow us to place constraints on the PAH abundance if the PAH band 
strengths are known. On the other hand, one can also infer information on the PAH band strengths from the 
aromatic absorption bands given knowledge of the PAH abundance derived from other considerations such 
as the 2175A hump and/or the 3 12/im PAH emission bands. For example, for fec=60ppm in PAHs (the 
maximum allowed by the 2175A hump), and assuming « 50 A, enhancement of the 6.2/im band by a factor 
Eq,2 > 3 would result in a 6.2/nm absorption feature which would be stronger than observed; conversely, if 
we had not enhanced the 6.2/um feature strength by a factor £^6.2=3, then the 6.2/im absorption feature 
strength in the model would be weaker than observed. However, this depends on the poorly-constrained 
value of aj, which determines the radius where carbonaceous grains make the transition from being PAH- 
like to becoming graphite-like. It is clear that (/ ArrfA~^)6.2 would increase (decrease) as is increased 
(decreased) (e.g., the = lOOA model results in (J ArrfA~^)6.2/Ar9.7 w 4.2 cm~^). Therefore, until we 
have better understanding of the PAH-graphite transition, wc should not over-state the constraints of the 
aromatic absorption bands on the PAH band strengths. However, the fact that our model is consistent with 
the 3.3, 6.2/im absorption bands, the 2175A hump, and the 3-12/im emission indicates that the choice of 
a^=50A is reasonable. 

We must also stress that the strengths of these PAH absorption bands could vary with physical conditions 

due to changes in the PAH ion fraction (/)ion(a)- In regions with increased cbion, the 6.2/im absorption feature 
would be strengthened and the 3.3/im feature would be weakened (see Table 1). The fact that there is 
no significant variation between the Galactic center sources and the local sources for the 6.2/tm absorption 
strength relative to the silicate optical depth (Schutte ct al. 1998) suggests that these rcigions probably have 
similar radiation field-to-electron density ratios. This is also indicated by the P(7.7/im)/i-'(11.3/tm) data of 
Chan et al. (2001) (sec their Figure 2c). 

Significant variations in the relative abundances of PAHs to large grains are expected in different envi- 
ronments due to depletion (e.g., accretion of PAHs into grain mantles, agglomeration of PAHs and other very 
small grains to form composite grains), production (e.g. fragmentation of graphitic grains, hydrogenated 
amorphous carbon [Scott, Duley, & Pinho 1997] or organic refractory materials [Grecnberg et al. 2000]), or 
drift of large grains relative to small grains due to anisotropic starlight (Weingartner & Draine 2001c). Such 
local variations in PAH abundance are indicated by the large cloud-to-cloud variations in mid-IR emission 
(Verter et al. 2000) as well as the lower UV extinction in dense regions. 

Ultrasmall silicate grains are excluded as a major component since the IRTS and ISO spectroscopic 
observations of the diffuse ISM show no evidence of the Si-0 9.7/im emission fcatiirc which would be expected 
for transiently-heated tiny silicate dust (Mattila et al. 1996; Onaka et al. 1996). However, the presence of 
these particles cannot be ruled out since a weak 9.7/tm emission feature may be hidden by the strong PAH 
features at 8.6/im and 11.3/im. On the basis of the observed IR emission spectrum for the diffuse ISM, 
the observed ultraviolet extinction curve, and the 10/im silicate absorption profile, Li & Draine (2001) have 
obtained upper limits on the abundances of ultrasmall (a < 15 A) amorphous and crystalline silicate grains. 
It is found that, contrary to previous work (< 1% of solar silicon abundance [Desert et al. 1986]), as much 
as ^10% of interstellar Si could be in a < 15 A silicate grains without violating observational constraints. 
Not more than ~5% of the Si can be in crystalline silicates (of any size). 
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12. Summary 

We have modeled the infrared emission spectra from dust in the diffuse ISM of the Milky Way, including 
the prominent emission features at 3.3, 6.2, 7.7, 8.6 and 11.3/Ltm and the continuum emission from near-IR 
to submillimctcr. Wc have compared our model to the observed emission from high galactic latitude regions 
as well as from regions in the Galactic plane. 

The dust model adopted here consists of a mixture of amorphous silicate grains and carbonaceous 
grains, both having a wide size distribution ranging from large grains ~l/im in diameter, down to molecules 

containing tens of atoms. We assume that the carbonaceous grains have graphitic properties at radii a > 50A, 
and polycyclic aromatic hydrocarbon (PAH)-likc properties at the very small sizes (see §2.2) which account 
for the 3.3, 6.2, 7.7, 8.6 and 11.3/im emission features seen in a wide range of objects. 

Based on recent progress made in both laboratory and theoretical studies of PAH absorption spectra, 
and guided by the astronomical observations, we have constructed "astronomical" absorption cross sections 
for both neutral and ionized PAHs from the far ultraviolet to the infrared (see §2.3, Table 1, and Appendix 
A) including dependence on grain size (which affects the onset of the absorption edge at long wavelength 
[see §A.2]), H/C ratio, and charge state. Detailed modelling indicates that the feature strengths of the 7.7, 
8.6(Um bands need to be increased by a factor ~2, and a factor -^3 for the 6.2/im band (see §2.3 and compare 
Figures 10a and 14a). In addition, a small amount of continuum absorption needs to be added to the PAH 
absorption spectrum to account for the near-IR continuum emission at l-bfim underlying the 3.3/xm C-H 
stretching feature detected in various sources (see §2.4). We represent this continuum opacity by adding 
about 1% of the opacity of graphite. We have also slightly modified the far-infrared and submillimeter 
emissivity of "astronomical silicate" in order to explicitly demonstrate that the observed emission from the 
interstellar medium can be reproduced using reasonable properties for carbonaceous and silicate grains. 

We model the heat capacities of carbonaceous and silicate grains using vibrational models developed 
by DLOl. The thermal-discrete method developed in DLOl is used to calculate the energy (temperature) 
distribution functions for transiently heated very small grains (see §6). We find that for grains heated 
by starlight with xmmp = 1, the energy distribution function is quite broad even for grains as large as 
a = 200 A, resulting in stronger emission at short wavelengths than would be computed using the "equilibrium 
temperature" T. The PAH charging computation is performed using the technique developed by Weingartner 
& Draine (2001b) (see §8). 

We adopt the size distributions for large grains (see §7) derived from modeling the interstellar extinction 
curve (Weingartner & Draine 2001a). We find that if the very small carbonaceous grain population includes 
two log-normal size distributions characterized by aoi = 3.5A, ao2 = 30A, and widths a = 0.4 for the PAH 
size distributions (with about 60ppm of C [relative to H] contained in PAHs; 45ppm in the former, and 
15ppm in the latter), then the overall grain model is able to closely reproduce the observed emission of 
the diffuse ISM from near-IR to submillimeter (see §10) while at the same time reproducing the observed 
extinction curve. Our dust model predicts a pronounced absorption feature at 6.2/xm and a much weaker 
feature at 3.3/xm due to PAHs with optical depths consistent with observations (see §11). 

We will show elsewhere (Draine & Li 2000b) that the very same ultrasmall grain component (i.e. 
PAHs) invoked here to account for the 3-60/im vibrational emission will produce "rotational" emission which 

can quite naturally explain the observed 10~100GHz "anomalous" microwave emission from the interstellar 
medium, confirming the original estimates by Draine & Lazarian (1998a, b). 

While real interstellar grains are no doubt more complex than the simple model presented here, this 
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extension of the original graphite-silicate grain model (Mathis, Rumpl, & Nordsieck 1977; Draine & Lee 1984) 
to explicitly include a PAH population leads to a well-defined grain model which is in excellent agreement 
with available observational constraints, ranging from the observed extinction curve to the observed infrared 
emission. 

We thank L.J. Allamandola, R. Arendt, .LM. Grcenberg, and J.C. Weingartner for helpful discussions. 
We are grateful to the anonymous referee for helpful comments. We thank D.M. Hudgins and L.J. Allaman- 
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Schlcgcl for help in obtaining the DIRBE fluxes for regions in the Galactic plane. We thank R.H. Lupton for 
the availability of the SM plotting package. This research was supported in part by NASA grant NAG5-7030 
and NSF grants AST-9619429 and AST-9988126. 



A. Optical Properties of Polycyclic Aromatic Hydroccirbons 

A.l. The UV and Far-UV Range 

Donn (1968) originally proposed that PAHs could contribute significantly to the interstellar UV extinc- 
tion. Although a single PAH species has strong UV absorption features (UV Atlas 1966), a cosmic mixture 
of many individual molecules, radicals, and ions, with many overlapping absorption bands, may effectively 
produce an absorption continuum, except perhaps in the visible, where some narrow features may stand out 
and account for some of the diffuse interstellar bands, and in the 2000-2400A region, where a concentration 
of strong absorption features may blend together to produce the 2175A extinction feature. 

Efforts over the past two decades to measure the absorption spectrum of PAHs have been primarily for 
individual ionized/neutral PAH molecules in hopes of identifying carriers of the diffuse interstellar bands (see 
Salama et al. 1996 and references therein). To the best of our knowledge, the only high quality absorption 
spectra published for natural mixtures of PAHs from the visible to far-UV are those of Leger et al. (1989b) 
and Joblin et al. (1992). For 3.3 ^m~^ < < 10/um~^, we adopt their experimental spectra as a starting 
point. All these spectra have a strong UV feature around 2175A (but relatively broader than the interstellar 
feature). We use the method of Fitzpatrick & Massa (1990) to fit the averaged spectra of Leger et al. 
(1989b)^ and of Joblin et al. (1992)^° by a Drude profile plus a linear continuum and a nonlinear far- 
ultraviolet (A"'^ > 5.9 /im~^) rise. To be consistent with the observed extinction, we then adjust the peak 
position of the ncar-UV hump (due to the tt-tt* transition) to be at A~i^. = 4.6 iim.~^ and the band width 
FWHM to be 7A^i^. — 1.0 /zm^^, with the integrated cross section conserved. We thus obtain eq. (10), 



"Wc adopted the PAH mixtures obtained from a coal pitch extract evaporated at 380K (averaged over those of 381K and 
385K), 410K, 440K (averaged over 436K and 440K), and 476K. In these mixtures over 300 molecules were identified with 
moleculaj- mass ranging from ~ 150 — 350 amu. The measurements were performed only for < 7.7 ixm~^, beyond which 
extrapolations were made. 

'^'^These mixtures were from a coal pitch extract evaporated at 570K, 600K and 630K, with molecular masses ranging from 
170 to 550 amu. 

^^The carrier of the 2175A extinction hump remains unidentified more than 35 years after its first detection (Stecher 1965). 
It is commonly attributed to some sort of graphitic carbonaceous material (see Draine 1989a). A correlation between the 2175A 
hump and the IRAS 12/im emission (dominated by PAHs) was found by Boulanger, Prevot, & Gry (1994) in the Chamaeleon 
cloud, suggesting a common carrier. Arguments against PAHs as the 2175A hump carrier can be found in Mathis (1998) where 
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for 3.3 < < 5.9^111-^ 

For A""'^ > 10 /im^^. no experimental data for PAH mixtures are available. We adopt that of coronene 
C24H12, characterized by a strong increase with A""'^ out to a strong peak at w 700A (due to the a-a* 
transition). Other simpler aromatic molecules such as naphthalene (CioHg) (Koch, Otto, & Radler 1972), 
benzene (CeHe) and its derivatives (Robin 1975) display similar spectral behaviour in the range of 10-30eV. 
For lOy^m"^ < A~^ < lb iiva~^, the laboratory measurements of coronene (Leger et al. 1989b) are fitted by 
a Drude profile plus a linear continuum to obtain eq. (7). 

For A^^ > 17.25 /im"-'^, we adopt the absorption cross sections per C atom calculated from the refractive 
indices of graphite (DL84); see eq. (5). A smooth transition from the absorption spectrum of coronene to 
that of graphite is made by extrapolating the coronene data to A~^ « 17.25 /um~^, resulting in eq. (6). 

To smoothly join our eq. (7) for A~^ > 10/im~^ to eq. (9) for A~^ < 7.7/um~\ for 7.7 < A~^ < 10/im~^ 
we adopt the polynomial eq.(8), which provides a good fit to the experimental data for PAH mixtures. 



A.2. The Visible and Near-UV Range 

The visible and near UV absorption properties of PAHs have been extensively studied (Clar 1964; Salama 
et al. 1996 and references therein). However, there are no high quality absorption spectra for PAH mixtures 
available in the literature. Therefore we adopt that of circumanthracene (C40H16), the largest symmetric 
PAH with the visible and near UV spectrum measured to date. The laboratory absorption spectrum of 
circumanthracene (Clar 1964) can be approximated by 

[abs(A)/A^c]c,oH,e = 2-54 X 10-17-3.431(A/;,m) ^^2/^. (Al) 

We adopt the cutoff function of Desert et al. (1990), 



cutoff (A, Ac) = — arctan 

TT 



lo'(y-i)' 
y 



+ \, y = Ac/A, (A2) 



where the cutoff wavelength Ac - the visual absorption edge - is derived from the extensive discussion of 
Salama et al. (1996): 



Ac 



/Ltm 



1/ (3.804M-0-5 + 1.052) , for neutral PAHs; , , 

1/ (2.282M-0-5 + 0.889) , for PAH cations; ^ ' 



where M, the number of fused benzenoid rings, is approximately M k, QANq for Nq > 40, M « 0.37Vc 
for A^c < 40 (A'c is the number of carbon atoms). Upon ionization, the visual absorption edge Ac shifts to 
longer wavelength; Ac also shifts to longer wavelength as the PAH size increases. 

To smoothly join the visual-near UV spectrum with that of A~^ > 3.3 ^m~^ (§A.l), we multiply the 
absorption cross section (Al) of circumanthracene by a factor 1.36, to obtain eq. (11) for 1 < A^^ < 3.3/xm~^. 



it was argued that the hump peak is sensitive to the PAH shape (but also see Duley & Scahia 1998) and tiie hump strength 
is senstive to the states of ionization (e.g. Lee & Wdowiak [1993] found that the 2175A hump is weakened in cations). In the 
dust model presented here the both the 2175A extinction hum and the 12//m emision are mainly contributed by PAHs, but we 
note that PAHs are not the major contributor to the far-UV (A < I8OOA) extinction. This is consistent with the observation 
that the far-UV extinction does not correlate with either the 2175A hump (Greenberg & Chlewicki 1983) or the IRAS 12^m 
emission (Boulanger et al. 1994). 
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Table 7. Observed PAH band widths (FWHM, iim) in selected sources. 



Objects 


3.3^m 


6.2^m 


7.7/^m 8.6/^ni 11.3/^ni 11.9/im 


12.7/im 


Reported res. 


Method" 


Spectrum sources 


diffuse ISM 


0.20'' 


0.67 


0.91 


0.55 


0.63 




AA = 0.3/nm 


Drude 


Onaka et al. 1996 


vdB 133 (little UV flux) 




0.30 


0.78 


0.50 


0.60 


0.85 


A/AA = 40 


Drude 


Uchida et al. 1998 


M17 (PDR) 




0.30 


0.68 


0.42 


0.60 


0.80 


A/AA = 40 


Drude 


Cesarsky et al. 1996a 


NGC 7023 (reflection nebula) 


0.04'= 


0.28 


0.60 


0.42 


0.43 


0.80 


A/AA = 40 


Drude 


Cesarsky et al. 1996b 


NGC 7023 (reflection nebula) 


0.04= 


0.17 


0.72 


0.28 


0.22 


0.33 


A/AA = 500 


Eye 


Moutou et al. 1999b 


NGC 1333 (reflection nebula)'^ 


0.04^ 


0.35 


0.70 


0.43 


0.40'' 




A/AA = 40 


Lorentzian 


Uchida et al. 2000 


p Oph (molecular cloud) 




0.27 


0.70 


0.42 


0.43 


1.18 


A/AA = 40 


Lorentzian 


Boulanger et al. 1996 


SMC 




0.20 


0.53 


0.36 


0.54 




A/AA = 40 


Lorentzian 


Reach et al. 2000 


Diffuse ISM 




0.20 


0.71 


0.42 


0.27 




AA = O.l^m 


Cauchy 


Mattila et al. 1996 


NGC 891 (external galaxy) 




0.22 


0.62 


0.30 


0.30 




AA = O.l/im 


Cauchy 


Mattila et al. 1999 


28 normal galaxies (averaged) 




0.36 


0.77 


0.45 


0.42 




AA = O.lfim 


Drude 


Hclou et al. 2000 


Orion ionization front 


0.048 


0.24'' 


O.8O'' 


0.28'' 


0.14 


0.27 


AX = 0.09/im 


Gaussian 


Roche et al. 1989 


NGC 2023 (reflection nebula) 


0.04° 


0.20 


0.74 


0.38 


0.36 




A/AA = 500 


Drude 


Verstraete et al. 2001 


Red Rectangle (post-AGB) 


0.0358 


0.16 


0.65 


0.36 


0.21 0.30 


0.42 


A/AA ^ 1000 


Eye 


Tielens et al. 1999 


30 Dor (LMC)' 


0.036 


0.30 


0.60 


0.40 


0.20 




A/AA = 1500 


Drude 


Sturm et al. 2000 


M82 (starburst galaxy) 


0.04 


0.20 


0.71 


0.42 


0.19 




A/AA = 1500 


Drude 


Sturm et al. 2000 


NGC 253 (starburst galaxy) 


0.042 


0.18 


0.78 


0.45 


0.19 




A/AA = 1500 


Drude 


Sturm et al. 2000 


Cirrinus(Soyforl 2 galaxv) 


0.045 


0.16 


0.82 


0.4:! 


0.20 




A/AA = 1500 


Drude 


Sturm cl al. 2000 


Present Model 


(1.0 


0.2(1 


0.70 


0. iO 


0.20'' O.:io 


o.:!0'- 






llu,. «'urk 



^FWHM values are determined by fitting the PAH features either with a set of Drude profiles, or Cauchy, Gaussian, Lorentzian profiles, 
or simply estimated by reading off the source spectra (hereafter "Eye"); those marked with "Cauchy", "Gaussian", or "Lorentzian" are 
taken from the source references; those marked with "Drude" or "Eye" are obtained in this work. 

^Tanaka et al. 1996 (resolution AA = O.lSjLtm). 

" Moutou ct al. 1999a (resolution A/AA ^ 500). 

^Uchida et al. (2000) have shown that the FWHM of the 6.2, 7.7, 8.6jum bands remains constant over different nebular positions with 
Go > 200; however, the 7.7yLim band appears to be broader at low Go- 

""Joblin ct al. 1996a (resolution AA ^ O-OOS^tm). 

^Joblin et al. 1996b (resolution A/AA ^ 60). 

^Geballe et al. 1989 (resolution AA = 0.009/Am). 

^Bregman et al. 1989 (resolution AA = 0.12j[im). 

^Except the 3.3 and 11.3/j-m bands, the signal-to-noisc ratio is very low for the 30 Dor spectrum (Sturm ct al. 2000). 

■'In a high resolution (A/AA 1400) survey of eight sources, Tokunaga et al. (1991) found that the 3.3/xm feature in extended objects 
such as planetary nebulae and HII regions has a band width of ^ 0.042/jm; this was further confirmed by Roche et al. (1996) who found 
18 planetary nebulae all have a band width of ^ 0.04//m for the 3.3/i.m feature (A/AA — 500). 

^On average, the lO-15/.im spectra of 16 HII regions, YSOs, and evolved stars (with resolution A/AA = 500 — 1500) show a band 
width of f« 0.20, 0.29)[/m for the 11.3, 12.7/im band, respectively (see Hony et al. 2001). 
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A.3. The Near-IR, Mid-IR Range 

During the last few years, a great deal of progress has been made by experimentalists and theorists in 
the study of the infrared properties of both neutral and ionized PAHs (Allamandola et al. 1999a; Hudgins 
& Allamandola 1999a, 1999b; Langhoff 1996; and references therein). The infrared spectra of PAHs are 
characterized by the C-C and C-H stretching/bending vibrational features. Observationally, the feature 
positions and widths are approximately constant in a wide variety of objects. Table 7 summarizes the band 
widths of the major PAH features observed in some well-observed astronomical sources. We represent the 
PAH infrared absorption (actually emission) cross sections by a set of Drude profiles with peak wavelengths 
\j and FWHM ^jXj taken from the observations^^ (see Table 7), respectively. Although the laboratory 
spectra of various PAH species are not exactly at the same wavelengths as the astronomical observations, we 
note that large symmetric compact PAHs are in close agreement with interstellar observations (see Langhoff 
1996). 

Although weaker than the 11.3/xm resonance, other C-H out-of-plane (oop) bending modes^^, have also 
been detected in various regions (Witteborn et al. 1989; Hony et al. 2001). From the point of view of 
stability, compact, symmetric structures are favored for interstellar PAHs (van der Zwet & Allamandola 
1985; Omont 1986; Allamandola ct al. 1989; Jochims ct al. 1994). Thus one expects a much lower quartet- 
CH (13.6yLtm) abundance than solo-, duet-, and trio-CH for interstellar PAHs, as has been recently confirmed 
by ISO observations (Hony et al. 2001). For large symmetric compact PAHs, it is reasonable to assume equal 
single, double, and triple adjacent CH units (Stein & Brown 1991). We therefore include Drude profiles with 
central wavelengths \j = 11.9, 12.7/im and FWHM jjXj = 0.30, 0.30/im in our PAH model spectrum for the 
duet-CH, trio-CH out-of-planc bending modes, respectively. 

Based on recent laboratory and theoretical studies of neutral and ionized PAHs (deFrees et al. 1993; 
Szczepanski & Vala 1993; Hudgins, Allamandola, & Sandford 1994; Hudgins & Allamandola 1995a, 1995b, 
1997; Langhoff 1996; Pauzat, Talbi, & EUinger 1997; Hudgins & Sandford 1998a, 1998b, 1998c) we estimate 
the integrated cross section of the C-C, C-H bands for neutral PAHs and PAH cations (normalized by 
the number of C atoms for C-C modes and by the number of H atoms for C-H modes). The results are 
presented in Table 1. We have considered only two charge states - neutral PAHs and PAH cations - since 
no experimental measurements have been made for multiply-ionized PAHs. It has been shown that the IR 
properties of PAH anions closely resemble those of PAH cations (e.g., see Szczepanski, Wehlburg, & Vala 
1995; Langhoff 1996; Hudgins et al. 2000) except for the very strong 3.3/zm C-H stretch enhancement in the 
anion^^. 



^^We rely on the highest resolution observations to indicate the FWHM of the narrower PAH features. 

i^The C-H oop bending modes are characteristic of the edge structure of PAHs. They occur at different wavelengths, 
depending on the number of neighbouring hydrogen atoms on one aromatic ring: for neutral PAHs, isolated C-H oop band 
(solo-CH) occurs at ~ 11.3/im, doubly adjacent CH (duct-CH) at ~ 11.9/im, triply adjacent CH (trio-CH) at ~ 12.7/im, 
quadruply adjacent CH at ~ 13.6/im. Upon ionization, the solo and duet CH modes shift to shorter wavelength (see Hudgins 
& Allamandola 1999b). 

^^Tlie .3.3;tm C-H stretching mode is about 3 times more intense in anions than the corresponding mode in neutral PAHs; 
however, no cation features have been observed in this mode (Hudgins et al. 2000). 
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A.4. The Far-IR Range 

The far-IR properties of PAHs are poorly known. Moutou et al. (1996) have carried out measurements 
of 40 PAH species in the range of 14 75/im. They found four prominent features at 16.2, 18.3, 21.2 and 
23.1/im (with band widths FWHM = 10, 20, 18, 20 cm^^, respectively) and a weak continuum beyond 
40/xm^^. To include these far-IR features, we adopt 4 Drude profiles for the 16.4, 18.3, 21.2 and 23.1/im 
bands and one Drude profile for the far-IR continuum (A > 24/im). Band widths are taken from Moutou et 
al. (1996) except for the 16.4/im feature, for which we adopt ~ 0.16/xm (Moutou et al. 2000). The integrated 
band strengths for the 16.4, 18.3, 21.2, and 23.1 /xm bands were estimated from the 40 species in Table 3 
of Moutou et al. (1996) by (1) summing the band strengths per C atom of features within ±20cm~^ of the 
nominal band centers, (2) dividing by 40, and (3) multiplying the resulting band strength per C atoms by 
a factor of 2 (a) since many of the species studied by Moutou et al. (e.g., hexabenzocoronene C42II18) arc 
highly symmetric and show only weak absorption in the 14-40/xm region and (b) in order to reproduce the 
observed strengthof the 16.4/im feature relative to shorter wavelength features (e.g., ll.Si^m) in the spectrum 
of NGC 7023 (Moutou et al. 1999a). 

The parameters for the FIR continuum Drude profile (<Tint,i4, A14, 714; see Table 1) were chosen by 
setting (Tint, 14 ~ Iii^m^c^^&hs^W'^^^^ ^ X]i=io^int.i averaged over the laboratory spectra of 40 PAH 
species (Moutou et al. 1996), and choosing A14, 714 to (1) minimize effects on the near-IR (A < 14/im) 
spectrum; (2) avoid a noticeable hump in the far infrared (A > 24/im) spectrum; (3) minimize absorption 
at A > Amax, where Amax corresponds to the lowest vibrational mode (see DLOl); (4) have the far infrared 
continuum dropping approximately as oc A~^ for A > 24/xm. 
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